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eo Abstract 
— In the last years it has been shown that some properties of strongly coupled 
Se superconductors can be potentially described by classical general relativity liv- 
eZ ing in one higher dimension, which is known as holographic superconductors. 
= This paper gives a quick and introductory overview of some holographic super- 
te conductor models with s-wave, p-wave and d-wave orders in the literature from 
"O : 3 . + . 
©. point of view of bottom-up, and summarizes some basic properties of these holo- 
ka graphic models in various regimes. The competition and coexistence of these 
PS superconductivity orders are also studied in these superconductor models. 
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1 Introduction 


The phenomenon of superconductivity was discovered in the early part of the last century 
that the electrical resistivity of a material suddenly drops to zero below a critical temper- 
ature Te. More importantly, the so-called Meissner effect tells us that the magnetic field 
is expelled in the superconducting phase, which is distinguished from the perfect conduc- 
tivity. In the latter case a pre-existing magnetic field will be trapped inside the sample. 
Conventional superconductors are well described by BCS theory [1], where the condensate 
is a Cooper pair of electrons bounded together by phonons. According to the symmetry 
of the spatial part of wave function of the Cooper pair, superconductors can be classified 
as the s-wave, p-wave, d-wave, f-wave superconductors, etc. However, some materials of 
significant theoretical and practical interest, such as high temperature cuprates and heavy 
fermion compounds, are beyond BCS theory. There are indications that the involving 
physics is in strongly coupled regime, so one needs a departure from the quasi-particle 
paradigm of Fermi liquid theory [2]. Condensed matter theories have very few tools to do 
this. 

On the other hand, although some of the deeper questions arising from the Anti-de 
Sitter/Conformal Field Theory(AdS/CFT) correspondence [8, 4, 5] remain to be under- 
stood from first principles, this duality creating an interface between gravitational theory 
and dynamics of quantum field theory provides an invaluable source of physical intuition 
as well as computational power. In particular, in a “large N and large A” limit, the grav- 
ity side can be well described by classical general relativity, while the dual field theory 
involves the dynamics with strong interaction.! It is often referred to as “holography” 
since a higher dimensional gravity system is described by a lower dimensional field theory 
without gravity, which is very reminiscent of an optical hologram. There are indeed many 
physical motivations that lead to this amazing holographic duality in the literature, such 
as renomalization group flow and black hole membrane paradigm. A very new perspec- 
tive was proposed, which is called the exact holographic mapping [6]. By constructing a 
unitary mapping from the Hilbert space of a lattice system in flat space (boundary) to 
that of another lattice system in one higher dimension (bulk), it provides a more explicit 
and complete understanding of the bulk theory for a given boundary theory and can be 
compared with AdS/CFT correspondence. 

It has been shown that the AdS/CFT correspondence can indeed provide solvable mod- 
els of strong coupling superconductivity, see refs. [7, 8, 9, 10] for reviews. The physical 
picture is that some gravity background would become unstable as one tunes some pa- 
rameter, such as temperature for black hole and chemical potential for AdS soliton, to 
developing some kind of hair. The emergency of the hair in the bulk corresponds to the 
condensation of a composite charged operator in the dual field theory. More precisely, 
the dual operator acquires a non-vanishing vacuum expectation value breaking the U(1) 
symmetry spontaneously. It has been uncovered that this simple holographic setup shows 


‘Loosely speaking, N? can be considered as the degrees of freedom in the dual field theory, and A as the 
characteristic strength of interactions. An elementary introduction to this correspondence can be found in 
the next section. 


similar properties with real superconductors. 

The holographic s-wave superconductor model known as Abelian-Higgs model was first 
realized in refs. [11, 12]. According to the AdS/CFT correspondence, in the gravity side, a 
Maxwell field and a charged scalar field are introduced to describe the U(1) symmetry and 
the scalar operator in the dual field theory, respectively. This holographic model undergoes 
a phase transition from black hole with no hair (normal phase/conductor phase) to the 
case with scalar hair at low temperatures (superconducting phase). The holographic model 
for the insulator/superconductor phase transition has been realized in ref. [13] at zero tem- 
perature by taking the AdS soliton as the gravitational background in the Abelian-Higgs 
model. Holographic d-wave model was constructed by introducing a charged massive spin 
two field propagating in the bulk [14, 15]. The superconductivity in the high temperature 
cuprates is well known to be of d-wave type. 

In recent years, evidence from several materials suggests that we now have examples 
of p-wave superconductivity, providing us new insights into the understanding of uncon- 
ventional superconductivity in strongly correlated electron systems [16]. To realize a holo- 
graphic p-wave model, one needs to introduce a charged vector field in the bulk as a vector 
order parameter. Ref. [17] presented a holographic p-wave model by introducing an SU(2) 
Yang-Mills field into the bulk, where a gauge boson generated by one SU(2) generator is 
dual to the vector order parameter. The authors of refs. [18, 19] constructed a holographic 
p-wave model by adopting a complex vector field charged under a U(1) gauge field, which is 
dual to a strongly coupled system involving a charged vector operator with a global U(1) 
symmetry. An alternative holographic realization of p-wave superconductivity emerges 
from the condensation of a two-form field in the bulk [20, 21, 22]. 

The philosophy for holographic setups is that even though the underlying microscopic 
description of the theory with a gravity dual is quite likely to be different form that 
arising in materials of experimental interest, it may uncover some universal aspects of 
the strongly coupled dynamics and kinematics, thus would help the development of new 
theories of superconductivity. By mapping the quantum physics of strongly correlated 
many body systems to the classical dynamics of black hole physics in one higher dimension, 
the holographic approach provides explicit examples of theories without a quasi-particle 
picture in which computations are nevertheless feasible. 

The models studied in the literature can be roughly divided into two classes, i.e., 
the bottom-up and top-down models. In the former approach the holographic model is 
constructed phenomenologically by picking relevant bulk fields corresponding to the most 
important operators on the dual field theory and then writing down a natural bulk action 
considering general symmetries and other features of dual system. Thus the holographic 
description is necessarily effective and can be used to describe a wide class of dual theories 
instead of a definite single theory. In the top-down approach, the construction of a model 
is uniquely determined by a consistent truncation from string theory or supergravity. One 
can usually have a much better control over the dual field theory, nonetheless, the resulted 
models are much more complicated. 

This paper aims at providing a quick and introductory overview of those three kinds of 


holographic superconductor models from the point of view of bottom-up.” The organization 
of the paper is as follows. In the next section, we review basic elements of Ginzburg-Landau 
theory of superconductivity and holographic duality, as a warm up. A brief introduction to 
the Abelian-Higgs model is presented in section 3. In section 4, we first introduce the SU(2) 
Yang-Mill model, focusing on its condensate and conductivity, then study the Maxwell- 
vector model, paying more attention to the vector condensate induced by magnetic fields 
and its complete phase diagram in terms of temperature and chemical potential, and finally 
discuss the third p-wave model by introducing a two-form in the bulk. This model can 
exhibit a novel helical superconducting phase. Section 5 is devoted to holographic d-wave 
models. In the next two sections, we pay attention to the competition and coexistence 
among different orders, including different superconducting orders in section 6 as well 
as superconducting order and magnetic order in section 7. The conclusion and some 
discussions are included in section 8. 


2 Preliminary 


2.1 Ginzburg-Landau theory 


The microscopic origin of traditional superconductivity is well understood by BCS theory, 
which explained the superconducting current as a superfluid of pairs of electrons interacting 
through the exchange of phonons. However, at phenomenological level, the Ginzburg- 
Landau theory of superconductivity [30] had great success in explaining the macroscopic 
properties of superconductors. ° 

In this phenomenological theory, the free energy of a superconductor can be expressed in 
terms of a complex order parameter field, Y, which is directly related to the density of the 
superconducting component. By assuming smallness of |W] and smallness of its gradients, 
the free energy near the superconducting critical temperature Te has the following form, 


|B)? 
2/0 
where Fn is the free energy in the normal phase, A is the vector potential and B = V x Ais 
the magnetic field. m* and e* are effective mass and charge of condensate. If one considers 
the BCS theory, m* = 2m and e* = 2e with m and e the mass and charge of electrons 
forming Copper pairs. a and p are two phenomenological parameters, which behave as 


a = ao(T. — T) with ap and 6 two positive constants. Note that we work in the unites 
h = c= 1. It is obvious that the free energy is invariant under the following transformation, 


\(-iV — e* A)U/? + 


F =F, talt + Siu + | (1) 


2m* 


s ee. al 
Y> Ue, A> A+ = V2), (2) 


?Holographic superconductor models constructed in the top-down approach can be found, for example, 
in refs. [23, 24, 25, 26, 27, 28, 29]. 
3To have a wide description of the Ginzburg-Landau theory see ref. [31] and references therein. 


which is known as the U(1) gauge symmetry. 
Minimising the free energy with respect to the order parameter and the vector potential, 
one obtains the Ginzburg-Landau equations, 


=> 


aY + B)YPW + 5 (iV — eA) Y = 0, (3) 

> ie* er? m 

J, = -— (WVU — UV") — —Y*y A. (4) 
2m* m* 


The first equation determines the order parameter Y, and the second one provides the 
superconducting current J, which is the dissipation-less electrical current. The Ginzburg- 
Landau theory actually can be derived from the BCS microscopic theory. Thus, the elec- 
trons that contribute to superconductivity would form a superfluid and ||? indicates the 
fraction of electrons condensed into a superfluid. 

This phenomenological theory can give many useful information even in homogeneous 
case. Let us consider a homogeneous superconductor with no superconducting current, so 
the equation for |W| simplifies to 


aY + pY = 0. (5) 


Above the superconducting transition temperature, T > Te, one only gets a trivial solution 
W = 0, which corresponds to the normal state of the superconductor. Below the critical 
temperature, T < T,, apart from the trivial solution, there are a series of non-trivial 


solutions which read 
a ao(T. — T) 
w,|=,/-_= . 6 
Vol =4/ 3 \/ 3 (6) 


Furthermore, compared with Y = 0, those solutions have lower potential energy, thus are 
dominant. Note that there are infinite solutions giving the ground state of superconducting 
phase. However, the true ground state can only choose one solution from them. Therefore, 
the ground state will change under the U(1) transformation (2). In such case, we call 
that the U(1) symmetry is spontaneously broken. From (6) one sees that Y approaches 
zero as T gets closer to T, from below, which is a typical behaviour of a second order 
phase transition. One can find from (4) that in the homogeneous case one can neglect the 
contribution from the first term and thus the superconducting current is proportional to 
the vector potential, i.e., J, œ A. If one takes a time derivative on both sides, one will 
obtain E = —0,A x 3J.. This means that the electric fields accelerate superconducting 
electrons resulting in the infinite DC conductivity. If one takes the curl and combines with 
Maxwell ’s equations, one will find V?B x B indicating the decay of magnetic fields inside 
a superconductor, i.e., the Meissner effect. 

The Ginzburg-Landau equations predict two characteristic lengths in a superconductor. 
The first one is the coherence length € which is given by 


| 1 
c= 2m*|al" (7) 
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It is the characteristic exponent of the variations of the density of superconducting com- 
ponent. In the BCS theory € denotes the characteristic Cooper pair size. The other one is 
the penetration length À which reads 


m* m* B 
Sy) — S| 8 
Loe*?|Wo|? \ Loe? |a| 8) 


where Wo is the equilibrium value of the order parameter in the absence of electromagnetic 
fields. This length characterises the speed of exponential decay of the magnetic field at the 
surface of a superconductor. 

Note that from definitions (7) and (8) the temperature dependences near T, behave as 
fot) (9) 
Ae T-T). (10) 
Both diverge as T — T, from below with the critical exponent 1/2. Nevertheless, the ratio 
k = X/E known as the Ginzburg-Landau parameter is temperature independent. Type- 
I superconductors correspond to cases with 0 < K < 1/./2, and type-II superconductors 
correspond to cases with « > 1/./2. One of the most important findings from the Ginzburg- 
Landau theory was that in a type-II superconductor, strong enough magnetic fields can 
penetrate the superconductor by forming the hexagonal lattice of quantised tubes of flux, 
called the Abrikosov vortex lattice. 

Finally one point we would like to emphasize is that in the Ginzburg-Landau theory the 
U(1) symmetry is broken spontaneously in the superconducting phase transition. Actually, 
only the spontaneous symmetry breaking feature itself can lead to many fundamental phe- 
nomenological properties of superconductivity, without any precise detail of the breaking 
mechanism specified [32]. In this review, we will show how a similar effective approach 
constitutes the basis of superconductivity in terms of holographic description. 


2.2 Holographic duality 


The original conjecture proposed by Maldacena [3] was that type-IIB string theory on 
the product spacetime AdS; x S° should be equivalent to M = 4 SU(N) supersymmetric 
Yang-Mills theory on the 3+1 dimensional boundary. This super-Yang-Mills theory is a 
conformal field theory, so this duality is named AdS/CFT correspondence. Later, this con- 
jecture has been generalized to more general gravitational backgrounds and cases without 
supersymmetry and conformal symmetry.’ From a modern perspective, the correspondence 
is an equality between a quantum field theory (QFT) in d dimensional spacetime and a 


4The simple examples are Lifshitz symmetry [33] and Schrödinger symmetry [34, 35], while more generic 
cases are those with generalized Lifshitz invariance and hyperscaling violation [36, 37, 38, 39], and the asso- 
ciated Schrödinger cousins [40]. However, those take us outside the best understood AdS/CFT framework. 
We shall focus on the most well defined case involving the bulk geometry with the asymptotically AdS 
behaviour. 


(quantum) gravity theory in d+1 spacetime dimensions. This correspondence is also some- 
times called gauge/gravity duality, gauge/string duality or holographic correspondence (or 
duality). 

A remarkable usefulness of the correspondence comes from the fact that it is a strong- 
weak duality: when the quantum field theory is strongly coupled, the dual gravitational 
theory is in a weakly interacting regime and thus more mathematically tractable, and vice 
versa. So the holographic duality provides us a powerful toolkit for studying strongly 
interacting systems. The backbone of the correspondence was elaborated by the authors 
of refs. [4, 5]. For every gauge invariant operator O in the QFT, there is a corresponding 
dynamical field ® in the bulk gravitational theory. The partition function in gravity side 
is equal to the generating functional of the dual boundary field theory. More specifically, 
adding a source J for O in the QFT is equivalent to impose a boundary condition for the 
dual field ® at the boundary of the gravity manifold (say at z — 0), i.e., the field ® tends 
towards the value 6 — ġo = J at the boundary up to an overall power of z. The formula 
reads 

Zoulk|® > bo = J] = (exp («| vm d’z 10) ) : (11) 

QFT 

where go is the determinant of the background metric of dual field theory. If we want to 
study a strongly coupled field theory, we can translate it into a weakly coupled gravity 
system. In the semiclassical limit, the partition function is equal to the on-shell action of 
the bulk theory and thus one only needs to solve particular differential equations of motion. 
Therefore we can compute expectation values and correlation functions of the operator O 
in the (strongly coupled) QFT by differentiating the left side with respect to J = ġo. 

In order to get familiar with the calculation by holography, let us consider O as a scalar 
operator which is dual to the bulk scalar field also denoted as ®. The minimal bulk action 
for ® is given by 


1 1 
oo = pare -300 — ie . (12) 
For illustration the gravity background is fixed as the pure AdS4+ı in Poincare coordinates 


2 


L 
ds? =< (d2 — di? + dë dz), (13) 


together with a profile for the scalar field ® = (z,t,#). Here # are d—1 spatial coordinates, 
t is a timelike coordinate and z is the radial spatial coordinate. L is known as AdS radius 
and the conformal boundary of AdS is located at z + 0. Note that the geometry is invariant 
under the scaling transformation (z,t,Z) > (Az, At, AZ). Actually, the full isometry group 
of AdSa+ı is identical to the conformal group in d dimensional boundary spacetime. 

To calculate the on-shell action for ®, we need to solve the equation of motion derived 
from action (12). Working in Fourier space ®(z,t,Z) => ®(z,w,k) = ®(z)e~**"*?, we 
obtain ne 

O26(z) — “—*a.0(2) m + ms )O(z)=0, k= —w? +P. (14) 


Near the boundary z — 0, the above equation admits the general asymptotic solution, 


=>. 


®(z,w,k) ~ Alk) 244 Bik) 24, 230, (15) 


with A = d/2 + \/m?L? + d2/4.5 The relation between A and B is determined by the 
interior of AdS. Making Fourier transformation back into real space, we then obtain 


D(z, t, Z) ~ AG, Z) 27 9+ B, Z) z^, 2 —0, (16) 


We now try to identify which term can be considered as the source J of the dual operator 
O. It turns out that, as long as m? L? > —d?/4 + 1, the mode A is non-normalizable with 
respect to the inner product 


(®,, ə) = -i | dzdť„/—gg"(®%0:®> = 020,07), (17) 
Lt 


with 4; a constant-t slice. The B mode in this case is normalizable. We identify the 
coefficient A as the source term, i.e., 


I(t, Z) = do(t,#) = A(t, #) = lim z^ D(z, t, Z). (18) 
This means that the on-shell action, and thus the partition function in AdS is a functional of 
J(t, Z). We can now argue that the scaling dimension of O is A without further calculation. 
Consider the scale transformation (z,t,#) > A(z, t, T), the scalar ® under such operation 
transforms as ®(Az, At, AZ) = @(z,t, 7). So the source term J = A must transform as 
J(At, AZ) = A4-4 I(t, Z), and thus according to (11) © transforms as O = A~4O which 
suggests the dimension of O should be A. 

Calculating the on-shell bulk action in terms of the solution with asymptotic expan- 
sion (16), one will find that the on-shell action will diverge near the boundary z > 0. This 
divergence is interpreted as dual to UV divergences of the boundary field theory. Actually, 
the infrared (IR) physics of the bulk near the boundary corresponds to the ultraviolet 
(UV) physics of dual QFT, and vice versa. This is called UV/IR relation [42] and the 
radial direction z plays the role of energy scale in the dual boundary theory. Physical 
processes in the bulk occurring at different radial positions correspond to different field 
theory processed with energies which scale as E ~ 1/z. 

The divergence can be cured by adding local counter terms at the boundary, known 
as holographic renormalization [43, 44, 45]. For the present case the counter term to be 
introduced is 


Sa = —— diry =y ®?, (19) 
2L z—0 
5Note that AdS 11 spacetime is stable even when the mass squared m? of scalar field is negative provided 
m? L? > mz,,L? = —d?/4 [41]. The lower bound mẹ p = —d?/4L? is often called Breitenloner-Freedman 
(BF) bound. 


where y is the determinant of the induced metric at the boundary. So the renormalized 
on-shell action should be S$” = So + Se. We can then compute the expectation value by 
using the basic formula (11). That relation implies 


_,Zourlt]  65°"[T] 


= 20 
where we have taken the semiclassical limit Zug = e5". Straightforward calculation 
shows that oA —d 

(O) (t, 2) = —>— BU 2). (21) 


This is often summarised as saying that the “non-normalizable” mode A gives the source 
in the dual field theory, whereas the “normalizable” mode B encodes the response. 

In the real world, many important experimental processes such as transport and spec- 
troscopy involve small time dependent perturbations about equilibrium. Those phenomena 
can be described by linear response theory, in which the basic quantity is the retarded 
Green’s function. The retarded Green’s function is defined to linearly relate sources and 
corresponding expectation values. In frequency space, it can be written as 


5(O) (w, k) = GE (w, k) 6J(w, k). (22) 


Using above formula one can continue to compute GË (w, k) which is given by 


(23) 


However, there is something subtle we shall discuss here. Different from the case in Eu- 
clidean signature where the bulk solution can be uniquely determined by additional re- 
quirement of regularity in the IR, while in the real-time Lorentzian signature, we must 
choose an appropriate boundary condition in far IR region of the geometry. This ambi- 
guity reflects multitude of real-time Green’s functions (Feynman, retarded, advanced) in 
the QFT. Since the retarded Green’s function describes causal response of the system to 
a perturbation, we involve an in-going condition describing stuff falling into the IR, i.e., 
moving towards larger z as time passes. The advanced Green’s function corresponds to 
the choice of out-going condition enforced in the IR region. ° 

Let us briefly consider the case with —d?/4 < m?L? < —d?/4 +1 where the second 
restriction comes from the unitary bound. One can easily check that both terms in (16) 
are normalizable with the inner product (17). So either one can be considered as a source, 
and the other one as a response. These two ways to quantize a scalar field in the bulk by 
imposing Dirichlet or Neumann like boundary conditions correspond to two different dual 
field theories [49], respectively. In the standard quantization, the corresponding operator 


®An intrinsically real-time holographic prescription was first proposed by the authors of ref. [46] by 
essentially analytically continuing the Euclidean prescription. It has been justified by a holographic version 
of the Schwinger-Keldysh formalism [47, 48]. 
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has dimension A, while in the alternative quantization, the corresponding operator has 
dimension d — A. 7 

The above discussion only uses the near boundary expansion (16) and thus applies 
to generic asymptotically AdS geometries. It can also be applied to other fields such as 
components of the metric and Maxwell fields. To sum up, we first obtain a solution which 
satisfies appropriate boundary conditions, especially the condition in the deep IR. Then 
we compute the properly renormalized on-shell action, identify the source and response 
from the asymptotic behaviour of the solution near the boundary, and compute the Green’s 
function through linear response. In the next section we will use this procedure to compute 
the optical conductivity. 

Another essential entry in the holographic dictionary is that the thermodynamic data of 
the QFT is entirely encoded in the thermodynamics of the black hole in the dual geometry. 
QFT states with finite temperature are dual to black hole geometries, where the Hawking 
temperature of the black hole is identified with the temperature in the QFT. Turning on 
a chemical potential in this QFT corresponds to gravity with a conserved charge. The 
thermal entropy of QFT is identified as the area of black hole horizon and the free energy 
is related to the Euclidian on-shell bulk action. As space is limited, we only introduce 
essential issues which will be needed to discuss holographic superconductors. 

Before the end of this subsection, let us point out that the holographic duality can be 
used to understand some hard nuts in quantum gravity from dual field theory side. A 
typical example is the black hole information paradox. It was first suggested by Hawk- 
ing [52] that black holes destroy information which seemed to conflict with the unitarity 
postulate of quantum mechanics. The black hole information paradox can be resolved, 
at least to some extent, by holography, because it shows how a black hole can evolve in 
a manner consistent with quantum mechanics in some contexts, i.e., evolves in a unitary 
fashion [53, 54]. There are some excellent papers talking about aspects of holographic 
duality, see, for example, refs. [55, 56, 57, 58, 59, 60] for more details. 


3 Holographic S-wave Models 


3.1 The Abelian-Higgs model 


In this subsection, we begin with the Abelian-Higgs model [11] by introducing a complex 
scalar field Y, with mass m and charge q, into the (3 + 1) dimensional Einstein-Maxwell 
theory with a negative cosmological constant. The complete action can be written down 


TIn fact, it has been shown that even more general quantisations are possible, like double trace defor- 
mation, see, for example, refs. [50, 51]. 
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as 
S= ga | ev= (r+ - = : wl’ — |VY — igAv/? — mop) : (24) 

where 2k? = 167G with G being the Newtonian gravitational constant, R is the scalar 
curvature of spacetime, L is the AdS radius and Maxwell field strength Fav = V,Av—VvAy- 
If one rescales A, + A,,/q and Y + W/gq, then the matter part has an overall factor 1/q? in 
front of its Lagrangian, thus the back reaction of the matter fields on the metric becomes 
negligible when q is large. The limit q —> oo with qA, and qY fixed is called the probe 
limit. Here we will review the results obtained in this probe approximation, which can 
simplify the problem while retains most of the interesting physics. The study including 
the back reaction of matter fields can be found in ref. [12]. 

The background metric is the AdS-Schwarzschild black hole with planar horizon 
dr? Cie me 2 2 3 [p83 
a a +dy"), f(r) =r —ry/r*), (25) 
where we have set the AdS radius L to be unity. The conformal boundary is located at 
r — oo. The Hawking temperature of the black hole is determined by the horizon radius 
rp: T = 3r;,/4a. The solution describes a thermal state of dual field theory in (2+1)- 
dimensions with temperature T. In addition, it is clear that the AdS-Schwarzschild black 
hole is an exact solution of the action (24) when the matter sector is negligible. It will 
be seen shortly that when the temperature is lowered enough, the black hole solution will 
become unstable and a new stable black hole solution appears with nontrivial scalar field. 

To see the formation of scalar hair, we are interested in static, translationally invariant 
solutions, thus we consider the ansatz [11] 


ds? = —f(r)dt? + 


W=uv(r), A=d(r)dt. (26) 


The r component of Maxwell equations implies that the phase of ~ must be constant. 
Therefore, for convenience, one can take % to be real. This leads to the equations of 
motion ° 


1 2 2 
Wae WE pli, 
F P F T 
" 2 / 2p? 
e ee 
a 


As pointed out in ref. [70], the coupling of the scalar to the Maxwell field produces a 
negative effective mass for ~ (see the third term in the first equation). Since this term 


8The model is a s-wave one since the condensed field is a scalar field dual to a scalar operator in 
the field theory side. This model can be straightforwardly generalized to other spacetime dimensions. 
Holographic s-wave superconductors with generalised couplings have also been considered in a number of 
works [61, 62, 63, 64, 65, 66, 67, 68, 69]. 

°In the probe limit, the concrete value of the charge q does not play an essential role. Without loss of 
generality, we take q to be one. 
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becomes more important at low temperatures, we expect an instability towards forming 
nontrivial scalar hair. 
The asymptotic behaviours of scalar field and gauge field near the AdS boundary are 


pote eye, pmp ten, (28) 


where Ay = (3 + V3? + 4m?)/2, u is the chemical potential and p is the charge density in 
the dual field theory. According to the AdS/CFT dictionary, the leading coefficient ~_ is 
regarded as the source of the dual scalar operator O with scaling dimension A,. Since we 
want the U(1) symmetry to be broken spontaneously, we should turn off the source, i.e., 
w_ = 0. Therefore the subleading term Y+ provides the vacuum expectation value (O) in 
the absence of any source. 1° 


T 
T 


Figure 1: The condensate as a function of temperature. The critical temperature T, is 
proportional to \/p. We choose m? = —2. 


Figure 1 shows how the condensate (O) behaves as a function of temperature in a 
canonical ensemble with p fixed to be one. As one can see that there is a critical temperature 
T, below which the condensate appears, then rises quickly as the system is cooled and 
finally goes to a constant for sufficiently low temperatures. This behaviour is qualitatively 
similar to that obtained in BCS theory and observed in many materials. Near the critical 
temperature T., (O) ~ (T, — T)'/?, which is the typical result predicated by Ginzburg- 
Landau theory, see equation (6). By comparing the free energy of these hairy configurations 
to the solution Y% = 0, ¢ = p(1/r;, —1/r) with no scalar hair, one finds that the hairy phase 
is thermodynamically favoured and the difference of free energies behaves like (T, — T)? 
near the critical point, indicating a second order phase transition. 

We now compute the optical conductivity, i.e., the conductivity as a function of fre- 
quency w, which is related to the retarded current-current two-point function for the U(1) 


10We only consider the standard quantization here, regarding the leading coefficient as the source of dual 
operator. An alternative way inducing spontaneous symmetry breaking in holographic superconductors is 
to introduce double trace deformation [71]. 
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symmetry, o(w) = EGP (w, k = 0). According to the holographic duality, this can be ob- 
tained by calculating electromagnetic fluctuations in the bulk. By symmetry, it is sufficient 
to turn on the perturbation 6A = A,(r)e~‘dz, then the linearized equation of motion for 
A, is 
A" j F y j a — 2y’ 
"of OP J 
To obtain the real time correlation functions for the dual boundary theory, the holo- 
graphic description associates in-going and out-going boundary conditions at the black 
hole horizon to retarded and advanced boundary correlators respectively [46]. To con- 
sider causal behaviour, one should impose the in-going wave condition at the horizon: 
A, ~ f~#/3'n, Near the AdS boundary, the asymptotic behaviour of A, is given by 


ASD (29) 


AW 
AgNO ee (30) 


According to the AdS/CFT correspondence, A) is the source, while A® is dual to the 
current. Thus one can obtain 


olw) = =G? (wv) = =—. (31) 


T J T T T 
1 oN | 
104 a ` - 
1 \ \ J 
1 ` \ J 
1 ` \ J 
O54 1 k J 
1 % ` ] 
1 i ` ] 
nite) 0.0 ee eects ie cane aCe ee 
“Lee oN ai ae ] 
\ ff \ we 4 
1 | ` 1 
L \ 1 
—0.5 + v; 1 1 
\ v 
7 \ I 
~V 
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Figure 2: The optical conductivity as a function of frequence. The solid lines in the 
left plot are the real part of the conductivity, while the dashed lines in the right plot 
are the imaginary part of the conductivity. We choose m? = —2. The horizontal lines 
correspond to temperatures above T}. Other curves from the left to right correspond to 
T/T. ~ 0.888 (blue), T/T, ~ 0.222 (green) and T/T. ~ 0.105 (red), respectively. There is 
a delta function at the origin for the real part of the conductivity in the condensed phase. 


The AC conductivity as a function of frequency is presented in figure 2. Above the 
critical temperature, the conductivity is a constant. As the temperature is lowered below 


Te, the optical conductivity develops a gap at some special frequency wg known as gap 
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frequency. As suggested in ref. [72], it can be identified with the one at the minimum of 
the imaginary part of the AC conductivity. Re[a(w)] is very small in the infrared and rises 
quickly at w,.'’ There also exists a small “bump” slightly above wg, which is reminiscent 
of the behaviour due to fermionic pairing [17]. For different choice of parameters, one 
can obtain a robust feature wg ~ 8T, with deviations of less than 10%. Compared to the 
corresponding BCS value w, ~ 3.57%, the result shown here is consistent with the fact that 
the holographic model describes a system at strong coupling. There is also a delta function 
at w = 0 appearing as soon as T < T}. This can be seen from the imaginary part of the 
conductivity. According to the Kramers-Kronig relation 


Im[o(w)] = -1p l ` PEs, (32) 


/ 
T o W =w 


one can conclude that the real part of the conductivity contains a Dirac delta function at 
w = 0 if and only if the imaginary part has a pole, i.e., Im(a) ~ 1/w. 

From above discussion, we see that this simple model can provide a holographically 
dual description of a superconductor. It predicts that a charged condensate emerges below 
a critical temperature via a second order transition, that the DC conductivity becomes in- 
finite, and that the optical conductivity develops a gap at low frequency. The temperature 
dependences of the coherence length € as well as the penetration length A in the holographic 
model are both proportional to (T, — T)~/? near the critical temperature [72, 74]. It has 
been shown that this holographic superconductor is type-II [12]. The condensate can form 
a lattice of vortices and the minimum of the free energy at long wavelength corresponds to 
a triangular array [75]. The effects of a superconducting condensate on holographic Fermi 
surfaces have been studied [76, 77]. All these features are very reminiscent of real supercon- 
ductors. Although the holographic model is very simple, it indeed captures some significant 
characteristics for superconductivity, thus helping us to understand real, strongly coupled 
superconductors. 


3.2 Holographic insulator/superconductor phase transition 


In this subsection, let us consider a five-dimensional Einstein-Abelian-Higgs theory with 
following action 


1 12 i ; 

S= T pew (r Po gret" — |V Y — iqA Y|? — m0?) . (83) 
When one does not include the matter sector, the theory has a five-dimensional AdS- 
Schwarzschild black hole solution. It is interesting to note that there also exists another 


"Tt has been shown that the conductivity is directly related to the reflection coefficient with the frequency 
given the incident energy [73]. The key point is that even as T — 0 there is still tunneling through the 
barrier provided by the effective potential. Therefore, a nonzero conductivity at small frequencies will 
always exist, and hence there is no hard gap in the optical conductivity at zero temperature. To obtain a 
superconductor with a hard gap, one might consider non-minimally coupled scalars in the bulk. 
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exact solution, so-called AdS soliton, in the theory (33). The AdS soliton solution can be 
obtained by double Wick rotation from the AdS-Schwarzschild black hole as 


2 
ds? = f(r)dy? + a + r?(—dt? + dr? + dy”), f(r) =r?(1—15/r*). (34) 
To remove the potential conical singularity, the spatial coordinate x has to be periodic 
with a period z/ro. If one considers the coordinates (x,r), the geometry looks like a cigar 
and the tip is given by r = ro. The AdS soliton has no horizon, and therefore no entropy is 
associated with this solution. Due to the existence of an IR cutoff at r = ro for the soliton 
solution, the field theory dual to this gravity background turns out to be in confined phase 
at zero temperature. Furthermore, this solution can be explained as a gravity dual to an 
insulator in condensed matter theory. If one increases the chemical potential to a critical 
value, the AdS soliton solution becomes unstable to developing a scalar hair with nontrivial 
scalar profile. It is shown that the new solution can describe a superconducting phase [13]. 
In this way, the holographic insultor/superconductor phase transition at zero temperature 
can be realized in the Abelian-Higgs model (33). 
More precisely, let us also consider the following ansatz in the probe limit 


Y= 4(r), Ay = o(r)dt. (35) 


In the AdS soliton (34) background, the equations of motions turn out to be 


" i m? "OA = 
w(t 7) -(F—‘ar)e=° 


(36) 
f J g Ww 
ots (E42 g -T4 =0. 
T f 
In the five-dimensional case, the BF bound is m% p = —4. For simplicity, let us consider the 
case with m? = —15/4. To solve the equations of motion, we have to specify the boundary 


conditions both at the tip and the AdS boundary. Near the AdS boundary, we have the 
following asymptotical form 


Note that in this case, both terms proportional to w_ and w are normalizable, so the 
corresponding operators QO, and O; have dimensions A = 3/2 and A = 5/2, respectively. 
On the other hand, near the tip of the soliton, these fields behave like 


w=atblog(r—ro)+c(r—1ro) +=, 


@=A+ Blog(r — ro) +CO(r— ro) +, (38) 


where a,b,c and A, B,C are all constants. The field regularity at the tip requires us 
to take b = B = 0. As in the previous subsection we can set q = 1 and further set 
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Figure 3: The behaviour of condensation for the operator (O1) (left) and (O2) (right) 
with respect to chemical potential. Used with permission from ref. [13]. 


ro = 1 without loss of generality. With the boundary conditions, solving the equations 
of motion, one can find that when the chemical potential u is beyond some critical value, 
the condensation happens. Concretely, for the operator O4, the critical chemical potential 
is uı = 0.84, while the critical chemical potential u> = 1.88 for the operator O2. The 
behaviour of condensation is plotted in figure 3 with respect to chemical potential. In 
figure 4 the charge density p with respect to chemical potential is plotted. We can see 
that at the phase transition point, its derivative is discontinuous, which verifies that the 
phase transition is indeed second order, since one has p = 02/0, where Q is the Gibbs 
free energy density. 


I I I 1 1 L a a a NC EE AE 
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Figure 4: The charge density p as a function of u when (O1) Æ 0 (left) and (Oz) 4 0 
(right). Its derivative jumps at the phase transition point. Used with permission from 
ref. [13]. 


To calculate conductivity o(w) we can consider the perturbation of the component A, 


in the soliton background. Assuming it has the form A, ~ e~’, its equation then turns 
out to be pa 2 opp? 
1 1 W q 
A —+-]A — — —— ] A, = 0. 39 
+ (547) et (eA) a 


At the tip one takes the Newmann boundary condition as in (38). Near the AdS boundary, 
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one has the asymptotical form as 


AP Aa? log Ar 
+ + 


A, = A® m 40 
where A is a cutoff. The holographic conductivity can be obtained as 
2AP iw 
olw) = -— + —. (41) 
wA® 2 


Since the background has no horizon, the real part of the conductivity always vanishes. 
This means that there is no dissipation. The imaginary part is plotted in figure 5: The left 
plot corresponds to the case of pure AdS solution without scalar hair, while the right one 
to the case with nontrivial scalar hair. There exist poles periodically at the points where 
A®) vanishes. These correspond to normalized modes dual to vector operators. One can 
see that when w is large, both case are similar, while when w —> 0, they are quite different. 
In the case without condensation, the imaginary part goes to zero when w —> 0, while it 
diverges in the case with condensation. According to the Kramers-Kronig relation (32), it 
shows that there is a delta functional support for the real part of conductivity at w = 0. 
Therefore the AdS soliton background with nontrivial scalar hair should be identified with 
the superconductivity. 


L bai a CO SO A ae Ce | 
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Figure 5: The imaginary part of the conductivity for the AdS soliton without condensation 
(O12) = 0 (left) and with condensation (O,) 4 0 (right). p = 0.0094 and u = 0.84 are 
taken in the right plot. Used with permission from ref. [13]. 


In the Einstein-Abelian-Higgs theory (33), besides the two phases described above, as 
in the four dimensional case, there exist another two solutions: AdS Reissner Nordström 
(AdS RN) black hole!” without scalar hair and AdS RN black hole with scalar hair, the 
latter can be identified with a superconductivity phase, while the former is dual to a con- 
ductor phase. Combining the four phases together, one could have the phase digram of 
the theory, which is schematically plotted in figure 6. The green line in the figure denotes 


12Tts precise form can be found in (65) below. 
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the first order phase transition, while two red lines represent second order phase transi- 
tion. Considering back reaction of matter sector, the complete phase diagrams in terms 
of temperature and chemical potential for the Abelian-Higgs model have been constructed 
in ref. [78]. It is interesting to note that the behaviour of the entanglement entropy with 
respect to chemical potential is non-monotonic and seems to be universal in this kind of 
insulator/superconductor models [79, 80, 81]. 


T 


AdS RN Black Hole 
1/n 


AdS BH 
Superconductor 


Figure 6: The phase diagram of AdS soliton and AdS black hole with a charged scalar 
field obtained in the large q limit. Used with permission from ref. [13]. 


4 Holographic P-wave Models 


4.1 The SU(2) Yang-Mills P-wave model 


The first holographic p-wave model is constructed by introducing a SU(2) Yang-Mills field 
in asymptotically AdS spacetime. One of three U(1) subgroups is regarded as the gauge 
group of electromagnetism and the off-diagonal gauge bosons which are charged under this 
U(1) gauge field are supposed to condense outside the horizon. The full action is given 
by [17] 

J= I diay gl- (R + $) - — Fe Fm] (42) 

2K? L? 402 ”” : 
where & is the four dimensional gravitational constant, g is the Yang-Mills coupling constant 
and L is the AdS radius. The field strength for the SU(2) gauge field A = Afir*dz* is 
Fin = OA, = OA, + er A AS, (43) 

where u,v = (t,r, x,y) denote the indices of spacetime and a, b,c = (1, 2,3) are the indices 
of the SU(2) group generators T° = o%/2i (o are Pauli matrices). e% is the totally 
antisymmetric tensor with e!?3 = +1. 
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Note that the ratio «/g measures the influence of Yang-Mills field on the background 
geometry. For the case «/g < 1 with A; fixed, the back reaction of the matter field can 
be ignored, thus the metric is simply AdS4 Schwarzschild black hole 


2 3 2 2 
ds? = a |- (1 — +) dt? + da? 4 iy? a (44) 
T 


with the temperature given by T = Hh, ‘ 


and we also fix a scale by setting ra = 1. 


Without loss of generality, we shall choose L = 1, 


4.1.1 Vector condensate 


To realize the p-wave condensate, one takes the ansatz 
A= O(r)r*dt + w(r)r'dz . (45) 


It is clear that the non-trivial profile of w(r) picks out the x direction as special, thus 
the condensed phase breaks the gauge group U(1)3 generated by 7? and SO(2) rotational 
symmetry in x — y plane. The relevant equations are [17] 


2 1 
o" += 6 — ap = 0, 
r r(r3 — 1) m 
M l + 2r? on r? o 0 l ) 
1 W ’ 
r(re—1)” ' @—12° © 
with primes representing the derivative with respect to r. 
The regularity at the horizon r = 1 demands the behaviour like 
d=d(r—-l1)+---, wa=uwotur(r—1)?4+---, (47) 
while the asymptotical expansion near the boundary r — oo takes the form 
W: 
d=p-E te, B= Wot ett (48) 


According to the holographic dictionary, u is regarded as chemical potential and p is the 
total charged density, and Wo is the source of the dual operator J”. To spontaneously break 
the U(1) symmetry, we should impose Wọ = 0, then the coefficient W; gives the vacuum 
expectation value of J”. According to the two-fluid model, the total charge density p can 
be divided into two components p = Pn + ps, where pn is the normal component, while ps 
is the superconducting component. In the holographic setup, the normal charge density pn 
is proportional to the r? part of the electric field at the horizon, i.e., pn = ¢,. Therefore 
the superconducting charge density is ps = p — Pn. 

By numerically solving the equations (46), one finds that the condensate is non-vanishing 
only when the rescaled temperature T’/,/p is small enough, i.e., lower than T, at which the 
condensate first turns on. As one can see in the right plot of figure 7, as the temperature 
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Figure 7: The fraction p,/p (right) and the condensate (left) as a function of temperature. 
Adapted with permission from ref. [17]. 


is lowered, (J) increases continuously. Near Te, (J”) vanishes as 7, — T, which is the 
typical behaviour predicted by Ginzburg-Landau theory. The fraction p,/p of the charge 
carried by the superconducting condensate goes to zero linearly near T,. 13 

We have interpreted U(1)3 generated by 7? as the gauge group of electromagnetism. 
The condensate of (J”} spontaneously breaks this U(1) symmetry as well as the rotational 
symmetry, thus resulting in an anisotropic superconducting phase. To see this much more 
clearly, we shall calculate the optical conductivity, which can be deduced by the retarded 
Green’s function of the U(1)3 current. Similar to the previous section, in gravity side the 
linear response to electromagnetic probes is turned out to study how linear perturbations 
of the T? component of the gauge field propagate. 


4.1.2 Conductivity 


In the presence of the condensate wr'dz, the x direction is special, so the conductivity oz, 
along the x direction is expected to be different from o,, along the y direction. To obtain 
consistent linearized equations, we can turn on the perturbation [17] 


A =e ™ [(apr! + afr*)dt + abr eda + aT” dy] ; (49) 


where all the functions depend on r only. By plugging the perturbation (49) into the 
linearized Yang-Mills equation, one finally obtains four second order equations 


2r3 +1 wr? T? 
3” 3/ 3 
y Tor Ln 21 rr = ay = 9, (50) 
r(r 1) (r 1) cr 1) 


13The ration ps/p versus temperature in the left plot is reminiscent of the temperature dependence of 
the superfluid of liquid He IT as measured from in the torsional oscillation disk stack experiment. However, 
we find ps goes to zero linearly here, while the experiment gives a critical exponent about 0.67. 
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aa AE a t gal -omal ial) 0 G) 
2 ow 
1⁄7 17 3 
| = 0, 51 
a, +7 + Ten (51b) 
1 2 y W : 
2" + =a = r(r3 = 1) (wa; + iwa’) = 0, (51c) 


and two first order constraint equations 


iwa? + daz — paz =0, 
1 1 (52) 
of 1 3) od 13 
iwa; — ġa; — (1 — 7g) 04, +¢a,+(1- 73)7 az =, 

It is clear that the equation of motion of the a mode decouples from the others, and the 
conductivity oy, exhibits similar “soft gap” behaviour to the s-wave model [11]. What 
we are interested in is the conductivity Css in the x direction. The conductivity oz, can 
be determined by solving the coupled equations (51) with the constraints given by (52). 
More precisely, we impose the ingoing wave condition at the horizon, which corresponds 
to a retarded Green’s function, 


az 


=(ro1y ee [1 + eM e140 = 17 +s -], 
al = (r — er Jae (r —1)7+ a, (r — 1)? +- | (53) 
a} = (r= 1) H/T Lat — 1) 4 O(n — 1)? ++], 


where all the coefficients can be fixed once wo, @, and w are specified. Near the conformal 
boundary r — oo, one has a generic solution to the equations of motion 


3(1) 
Wag gist, 
(0) , At” a0) _, Ar” ea 
(SAS) p p SA or a a 


As pointed out in ref. [17], there exists a residual gauge invariance. After fixing this residual 
gauge freedom, one can finally obtain the gauge invariant conductivity along x direction 


wA 4 a 


(55) 


pe? — w? 


-Í [gu 
Fae AO) [Ls a 


l4However, by considering the back reaction to the metric in the SU(2) model (42), it has been shown 
that the conductivity in the y direction has a “hard gap” at zero temperature, i.e., the real part of the 
conductivity is zero for an excitation frequency less than the gap frequency [82]. 
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Numerical calculation can only display the continuous part of o,.(w). One can reveal the 
non-analytic behaviour by virtue of the Kramers-Krong relations, which tells us that a sim- 
ple pole in Im[o,,(w)] at wo implies a delta function 8(w — wo) to Re[o,.(w)]. Further more, 
the positivity constraint on the real part of conductivities requires any pole of Im[o,.(w)| 
on the real axis to have a positive residue. 

The behaviour of conductivities as a function of frequency w is shown in figure 8, from 
which one can see the following features [17]. First, both czs and oyy approach constant for 
sufficiently large w. This is because the condensate involves dynamics with a characteristic 
energy scale set by ,/p. If w > ,/p, the propagation of the gauge boson should become 
insensitive to the condensate and can be approximated by the case in pure AdS4, thus is 
a constant. Second, o,, exhibits gapped dependence similar to the Abelian-Higgs model 
in figure 2. Re(o) is very small in the infrared, then rises quickly at w = wg ~ \/p. There 
is a slight “bump” a little above w, which is reminiscent of the behaviour expected for 
fermionic pairing. Third, there is a pole in Im[oz,] at w = wo ~ 1.8,/p. Therefore, there is 
a delta function contribution to Re[oz,,| at w = wo. Finally, in the small w region, Re[oz,| 
can be well parameterized in terms of the Drude model 


00 
1+ wr?’ 


Rel[oprude] = (56) 
where 09 gives the DC conductivity and 7 is the scattering time. The best fit gives a 
narrow Drude peak in o,, and suggests conductivity due to quasi-particles with scattering 
time to diverge as T > 0. 
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Figure 8: Conductivities oz, and oy, with respect to frequency at T’/,/p œ 0.0779. The 
dotted curves are the best fits of the Drude model prediction to Re o,,(w). Adapted with 
permission from ref. [17]. 


We do not have a microscopic description of the condensate in the language of the dual 
theory without gravity. However, we know clearly that there is an SU(2) current algebra, 
and the component J” develops an expectation for sufficiently large chemical potential. 
Yet we only turn on t'dx mode corresponding to the p-wave background. The (p + ip)- 
wave case can be realized by involving the combination t'dx + r?dy. This mode results 
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in an isotropic superconducting phase which exhibits a pseudogap!” at low temperatures 
and a nonzero Hall conductivity with no external magnetic field [83]. However, it should 
be pointed out that (p + ip) configurations are unstable against turning into pure p-wave 
background. The insulator/superconductor phase transition for the SU(2) p-wave model 
has been studied in ref. [84]. 

A new ground state can be found when a magnetic component of the gauge field is larger 
than a critical value, which forms a triangular Abrikosov lattice in the spatial directions 
perpendicular to the magnetic field [85, 86]. In the same spirit, a p-wave superconductor 
for which the dual field is explicitly known has been constructed in refs. [87, 88, 89] by 
embedding a probe of two coincident D7-branes in the AdS black hole background. From 
this top-down approach one can try to identify the SU(2) chemical potential as an isospin 
chemical potential and the condensate as a p meson. The back reaction of the gauge field on 
the metric in the SU(2) Yang-Mills model has been considered in refs. [90]. It is interesting 
to note that when the back reaction is strong enough, the phase transition will be a first 
order one. The holographic SU(2) p-wave superconductor model has been extended to 
include, for example, the Gauss-Bonnet term [91, 92] and Chern-Simons coupling [93]. In 
addition, based on the backreacted metric, the behaviour of entanglement entropy in the 
holographic superconducting phase transitions has been studied in refs. [94, 95, 96]. 


4.2 The Maxwell-Vector P-wave model 


Let us introduce a charged vector field into the (d+1) dimensional Einstein- Maxwell theory 
with a negative cosmological constant. The full action reads [18, 19] 


1 


d(d— 1 
=r | da/—-g(R+ Ac) + £m); 
2K L (57) 
1 ar oe 7 
bm = — i — aad =m pho" + igypupLF™, 


where a dagger denotes complex conjugation and p, is a complex vector field with mass m 
and charge q. We define ppv = D,p, — Dup, with the covariant derivative D, = Vp — iqAp. 
The last non-minimal coupling term characterizes the magnetic moment of the vector field 
Pu: 

Since p, is charged under the U(1) gauge field, according to AdS/CFT correspondence, 
its dual operator J” will carry the same charge under this symmetry and a vacuum expec- 
tation value of this operator will then trigger the U(1) symmetry breaking spontaneously. 
Thus, the condensate of the dual vector operator will break the U(1) symmetry as well as 
the spatial rotational symmetry since the condensate will pick out one direction as special. 
Therefore, viewing this vector field as an order parameter, the holographic model can be 
used to mimic a p-wave superconductor (superfluid) phase transition. The gravity back- 
ground without vector hair (p, = 0)/with vector hair (p, 4 0) is used to mimic the normal 
phase/superconducting phase in the dual system. 


15The terminology “pseudogap” here is to denote a well defined gap in the dissipative conductivity at 
low frequencies in which the conductivity is not identically zero. 
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Indeed, it was shown in ref. [18] that working on the probe limit, as one lowers the 
temperature, the normal phase becomes unstable to developing nontrivial configuration of 
the vector field. The calculation of the optical conductivity reveals that there is a delta 
function at the origin for the real part of the conductivity, which means the condensed phase 
is indeed superconducting. In this subsection, we shall review the effect of a background 
magnetic field on the model and its complete phase diagram in terms of temperature and 
chemical potential. 


4.2.1 Condensate induced by magnetic field 


Generally speaking, to consider the case with a magnetic field, one needs to solve coupled 
partial differential equations which is much more involved in practice. However, if one is 
interested in the instability induced by the magnetic field, one can overcome this difficulty 
by only focusing the dynamics near the critical point at which the condensate is very small. 
More precisely, one can introduce a deviation parameter e€ from the critical point at which 
the condensate begins to appear. The coupled equations of motion can then be solved 
order by order in terms of the power of e. 

Following the above procedure, we now turn on a magnetic field to study how the 
applied magnetic field influences the system. The background is taken to be a (3+1) 
dimensional AdS-Schwarzschild black hole (25). A consistent ansatz is as follows [18] 


pdz” = [ep,(r, x)e™ + O(e?)]dx + [epy(r, zee” + O(e*)]dy, 


Avda” = [ġ(r) + O(®)]dt + [Ba + O(®)]dy, (58) 


where pz(r,2), py(r,x) are all real functions, p is a real constant and the constant 0 is 
the phase difference between the x and y components of the vector field p,. The constant 
magnetic field B is perpendicular to the x — y plane. 
The profile of ọ can be uniquely determined at the zeroth order of e, which takes the 
form 
olr) =e — r/r) , (59) 
with u interpreted as the chemical potential. The equations of motion for p, and p, can be 
deduced from (57) at order O(¢). We further separate the variables as p,(r, x) = yz(r)X(z) 
and p,(r, x) = Y,(r)Y (x). Then one can get the following equations '° 
1 i J MOA m2 E 
Pr Pr — Pr — Pu 
pe gee Ee ey 
—X F (1+~7)qBY + (qBx — p}? X = EX, 
~Y + (1+ 7)qBX + (qBx — p)’Y = EY, 


ph + =0, 


(60) 


where the prime denotes the derivative with respect to r and the dot denotes the derivative 
with respect to x. We have also made a consistent assumption Yy = Yy and E is a constant 


16 Tn order to satisfy the equations of motion with the given ansatz, 0 can only be chosen as 0, = 3 +2n7 
or 0- = —5 + 2nz with n an arbitrary integer. Here and below the upper signs correspond to the #4 case 
and the lower to the 6_ case. 
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coming from variables separation. The last two equations for X (x) and Y (x) can be solved 
analytically and the eigenvalue is given by E = (2n + 1)|q4B| + (1 + 7)qB where n can be 
chosen as a non-negative integer. 

We are interested in how the applied magnetic field influences on the transition tem- 
perature from the normal phase to the condensed phase. The effective mass of the charged 
vector field in the lowest energy state, i.e., in the lowest Landau level n = 0 depends on 
the magnetic field B and the non-minimal coupling parameter y as 


2 2 |yqBl Ge 
in = mMm = —— = ——. 
eff r2 f 


It is clear that the increase of the magnetic field B decreases the effective mass and thus 
tends to raise the transition temperature, even in the case that the electric field is turned off. 
Only the magnetic field itself can trigger the phase transition. This result has an analogy to 
the QCD vacuum instability induced by a strong magnetic field to spontaneously developing 
the p-meson condensate. It is clear that the last term in (57) describing a non-minimal 
coupling of the vector field p“ to the gauge field A,, plays a crucial role in the instability. 
Note that similar coupling can be found in many formalisms used to describe the coupling 
of magnetic moment to the background magnetic field for charged vector particles [97, 98]. 
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Figure 9: The transition temperature from the normal phase to the condensed phase as 
a function of magnetic field. T, is the critical temperature in the case without magnetic 
field. The magnetic field raises the transition temperature. One has chosen m? = 3/4. 
This figure was taken from ref. [18]. 


The (T, B) phase diagram for the lowest Landau level is depicted in figure 9 in the case 
with fixed charge density p = urn. To determine which side of the phase transition line is 
the condensed phase, we can consider the equation (61). It suggests that the magnetic field 
decreases the effective mass. So if we increase the magnetic field at a fixed temperature, 
the normal state will become unstable for sufficiently large magnetic field. 

It is clear that the transition temperature increases with the applied magnetic field. 
In ordinary superconductors an external magnetic field suppresses superconductivity via 
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diamagnetic and Pauli pair breaking effects. However, it has also been proposed that 
the magnetic field induced superconductivity can also be realized in type-II supercon- 
ductors [99, 100], in which the Abrikosov flux lattice may enter a quantum limit of the 
low Landau level dominance with a spin-triplet pairing. And possible experimental ev- 
idence for the strong magnetic induced superconductivity can be found, for example, in 
refs. [101, 102). 

Due to the degeneracy in p, a linear superposition of the solutions with different p is 
also a solution of the model at O(e). We can take this advantage to construct a class of 
vortex lattice solutions. As a typical example, the triangular lattice is shown in figure 10. 
It should be stressed that it is the special combinations J+ = (J* + ifv} which exhibit 
the vortex lattice structure. Strictly speaking, to obtain the true ground state, one should 
calculate the free energy of the solutions with different lattice structures from the action 
to find which configuration minimizes the free energy. It turns out that the linear analysis 
presented here is not sufficient to determine the most stable solution, thus should include 
higher order contributions. Furthermore, it is worthwhile to mention that in the AdS 
soliton background, the external magnetic field triggered phase transition and vortex lattice 
structure also happen for the vector field p-wave model [103]. 


Figure 10: The vortex lattice structure for the triangular lattice in z — y plane. The 
contour plot is also drawn in the bottom. In particular, the condensate vanishes in the 
core of each vortex. The figure was taken from ref. [18]. 


The response of this system to the magnetic field is quite different from the behaviour of 
ordinary superconductor where the magnetic field makes the transition more difficult. But 
the result here is quite similar to the case of QCD vacuum instability induced by strong 
magnetic field to spontaneously developing the p-meson condensate [104, 105]. Although 
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so, it was shown that in model (57) the condensate of the vector operator forms a vortex 
lattice structure in the spatial directions perpendicular to the magnetic field. Of course, 
the non-minimal coupling term in the action plays a crucial role in both cases. Therefore 
in some sense, this model is a holographic setup of the study of p-meson condensate. 


4.2.2 The complete phase diagram 


The probe approximation neglecting the back reaction of the matter fields can indeed 
uncover many key properties. Nevertheless, it still loses some important information, such 
as the phase structure of the system. In the following paragraphs, we will discuss both the 
black hole background and soliton background in full back reaction case. Then a complete 
phase diagram in terms of temperature and chemical potential will be shown. We shall 
consider a (4 + 1)-dimensional bulk theory [106]. 

We would like to study a dual theory with finite chemical potential or charge density 
accompanied by a U(1) symmetry, so we turn on A; in the bulk. We want to allow for states 


A 


with a non-trivial current (J,,), for which we further introduce p, in the bulk. Because a 
non-vanishing (J,) picks out x direction as special, which obviously breaks the rotational 
symmetry in spatial plane, thus we should introduce an additional function in the xx 
component of the metric in order to describe the anisotropy. Therefore, for the matter 


part, we consider the ansatz 
pdz” = pz(r)dxr, A, dz” = ġ(r)dt. (62) 


We will consider black hole and soliton backgrounds separately. 


(1) AdS black hole with vector hair. 
For the black hole background, we adopt the following metric ansatz 


d 2 

ds? = —a(rje dt? + T + r° (c(r)dz? + dy? + d2?) . (63) 
a(r 

The position of horizon is denoted as r, at which a(r) = 0 and the conformal boundary is 
located at r — oo. One finds that the r component of Maxwell equations implies that the 
phase of p, must be constant. Without loss of generality, we can take p, to be real. Then, 
the independent equations of motion in terms of above ansatz are deduced as follows 
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where the prime denotes the derivative with respect to r. 
When py = 0, there exists an exactly analytical black hole solution, namely, AdS 
Reissner-Nordstrom black hole which reads 
d 2 
SS + 12(da? + dy? + dz), 
f(r) 
2 


2 we Tr\? H eR? TRN? 

a h- (1+ 5) (>) + 3p a f s) =n (*) |. 
This solution is dual to a conductor phase in the dual field theory. However, the full coupled 
equations of motion do not admit an analytical solution with non-trivial p,. Therefore, we 
have to solve them numerically. We will use shooting method to solve equations (64). In 
order to find the solutions for all the five functions, i.e., p.(r), ọ(r),a(r),b(r) and c(r) one 
must impose suitable boundary conditions both at conformal boundary r — oo and at the 
horizon r = r}. 


In order to match the asymptotical AdS boundary, the general falloff near the AdS 
boundary behaves as 


ds? = —f(r)dt? + 
(65) 


P Pa 
C= g tts Pane 


b 
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Pr+ 
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(66) 


where the dots stand for the higher order terms in the expansion of 1/r and A4 = 1 
JV1+m2.‘" In general, in the above expansion we must impose p,_ = 0, which meets 
the requirement that the condensate appears spontaneously. According to the AdS/CFT 
dictionary, up to a normalization, the coefficients jz, p, and pz, are regarded as chemical 
potential, charge density and the x component of the vacuum expectation value of the 
vector operator J” in the dual field theory, respectively. 

We focus on black hole configurations that have a regular event horizon located at 
r, and require the regularity conditions at the horizon r = ra, which means that all five 
functions {pz, ¢, a, b, c} would have finite values at r, and admit a series expansion in terms 
of (r— rn). After substituting such series expansion into equations (64), one finds there are 
only six independent parameters at the horizon, i.e., {rp, P2(Tn), 94ra), (rn), (rn) } and 
other coefficients can be expressed in terms of those parameters. 

Two free parameters b(r;,) and c(r,) can be fixed by AdS boundary conditions that 
b(r > co) = 0 and c(r > co) = 1. Without loss of generality, the location of r, can be 
fixed to be one in our numerical calculation. We are then left with two independent param- 
eters {px(rn), '(rn)}. By choosing ¢'(r;,) as the shooting parameter to match the source 
free condition, i.e., pr = 0, we finally have a one-parameter family of solutions labeled 
by the value of p, at the horizon. After solving the set of equations, we can read off the 


17The m? has a lower bound as m? = —1 with A, = A_ = 1. In that case, there exists a logarithmic 
term in the asymptotical expansion of py. One has to treat such a term as the source set to be zero to 
avoid the instability induced by this term [72]. We will always consider the case with m? > —1. 
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condensate (J z), chemical potential u and charge density p directly from the asymptotical 
expansion (66). 


(2) AdS soliton with vector hair. 
To construct homogeneous charged solutions with vector hair in the soliton background, 
we take the metric as 
d 2 
ds? = Fay + (Flr at? + Merde? + dy? + g(r) dy?) (67) 
r2g(r 
where g(r) vanishes at the tip r = ro of the soliton. The asymptotical AdS boundary is 
located at r — oo. Further, in order to obtain a smooth geometry at the tip ro, n should 


be made with an identification 


x(ro) 
4te 2 


- rêg'(ro) ` 
This gives a dual picture of the boundary theory with a mass gap, which is reminiscent of 


an insulating phase. 
The independent equations of motion are deduced as follows 


nvrntl, (68) 
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where the prime denotes the derivative with respect to r. Similar to the black hole case, 
we will solve those coupled equations of motion numerically by use of shooting method. In 
order to find the solutions for all the six functions F = {pz,¢, f, g, h, x} one must impose 
suitable boundary conditions at both conformal boundary r — oo and the tip r = rọ. 

The asymptotical expansion for metric fields and matter fields near the boundary r > 
oo is as follows 


B pP _ Pe- Pat 2 fa 
d a Pz = ke A pair at 
g X ey) 
4 4 4 
eas i h=1l+ at 3 a j 


where the dots stand for the higher order terms of 1/r. We choose the source free condition 
Px- = 0 as before. The coefficients u, p, and pr}, are directly related to the chemical 
potential, charge density and x component of the vacuum expectation value of the vector 
operator J” in the dual system, respectively. 

We impose the regularity conditions at the tip r = ro, which means that all functions 
have finite values and admit a series expansion in terms of (r — ro) as 


F = Fro) + F"(ro)(r — 10) +++ - (71) 


By plugging the expansion (71) into (69), one can find that there are six independent 
parameters at the tip {r0, Pz(ro), (ro), f (ro), A(ro), x(ro)}. However, there exist four useful 
scaling symmetries in the equations of motion, which read 


XOxX+A, >een, (72) 
ġ > id, t> AH, fr 7, (73) 
Pz > àp, >A r, h> Xh, (74) 
and 
r — Ar, (t,£, y, n) —> ATHE, £,Y, n), ($, Px) > AQ, Px) , (75) 


where in each case A is a real positive constant. 

By using above four scaling symmetries, we can first set {ro = 1, f(ro) = 1, h(ro) = 
1, x(ro) = 0} for performing numerics. After solving the coupled differential equations, one 
should use the first three symmetries again to satisfy the asymptotic conditions f(co) = 1, 
h(oo) = 1 and y(co) = 0. We choose (ro) as the shooting parameter to match the 
source free condition, i.e., ps_ = 0. Finally, for fixed m? and q, we have a one-parameter 
family of solutions labeled by p,(r9). After solving the set of equations, we can read off 
the condensate (Je), chemical potential u and charge density p from the corresponding 
coefficients in (70). It should be noticed that different solutions obtained in this way will 
have different periods I for 7 direction. We should use the last scaling symmetry to set all 
of the periods I equal in order to obtain same boundary geometry. We shall fix I to be m 
in this section. 

These two kinds of situations have been well studied in ref. [106]. There are four different 
bulk solutions given by the pure AdS soliton, AdS Reissner-Nordstrom black hole and their 
vector hairy counterparts. According to the AdS/CFT dictionary, the hairy solution is dual 
to a system with a non-zero vacuum expectation value of the charged vector operator which 
breaks the U(1) symmetry and the spatial rotation symmetry spontaneously. The above 
four solutions in the bulk correspond to an insulating phase, a conducting phase, a soliton 
superconducting phase and a black hole superconducting phase, respectively. Since we do 
not turn on magnetic field, the model is left with two independent parameters, i.e., the 
mass m of the vector field giving the scaling dimension of the dual vector operator and 
its charge q controlling the strength of the back reaction on the background geometry. 
The phase structure of the model heavily depends on those two parameters. There exist 
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Figure 11: The complete phase diagrams of the Maxwell-vector model with S=pure AdS 
soliton, BH=AdS Reissner-Nordström black hole, SC=hairy soliton, and BC=hairy black 
hole. In each region the thermodynamically stable phase is labeled. As m? and q are 
changed, the shape of each region gets modified. This figure is described in ref. [106] in 
more detail. 
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second order, first order and zeroth order! phase transitions as well as the “retrograde 
condensation” in which the hairy solutions exist only above a critical temperature or below 
a critical chemical potential with the free energy much larger than the solutions without 
hair. 

With four kinds of phases at hand, the complete phase diagrams can be constructed in 
terms of temperature and chemical potential. At each point in T-u plane, one should find 
the phase which has the lowest free energy. Since there are many types of phase transitions 
in both soliton and black hole backgrounds, the T-u phase diagrams are expected to be 
much more complicated than the holographic s-wave model [78] and the Yang-Mills p- 
wave model [84]. Some typical examples are shown in figure 11. We can see from the 
complete phase diagrams that in some cases, more than one superconducting phase appears 
in a phase diagram in the model. The phase diagrams for some realistic superconducting 
materials are usually complicated, and indeed, more than one superconducting phase can 
occur, for example, see refs. [111, 112, 113]. Definitely, it is of great interest to see whether 
this model is relevant to those superconducting materials. 


4.3 The Helical P-wave model 


The gravity solutions above mainly describe spatially homogeneous superconducting states. 
However, it has long been known that it is possible to have superconducting states that are 
spatially inhomogeneous. A well known example is the Fulde-Ferrell-Larkin-Ovchinnikov 
(FFLO) phase, for which a Cooper pair consisting of two fermions with different Fermi mo- 
menta condenses leading to an order parameter with non-vanishing total momentum [114, 
115]. In this section, we shall introduce a holographic model which can realize p-wave 
superconducting phase with a helical order. That is to say, the order parameter points in a 
given direction in a plane which then rotates as one moves along the direction orthogonal 
to the plane. 

We consider a (4+1) dimensional model with a gauge field A, and a charged two-form 
Oi (el; 22) 

j= f Px /—g9[R +12- Fwy FH = F 2. sagt Cw Ho] , (76) 

where we have chosen units where the AdS radius is unity, a dagger denotes complex 
conjugation and the field strengths read 


F=dA, H=dC+ieAAC. (77) 


The gauge field A, is dual to a current in the dual theory and the two-form Cy, 
corresponds to a self-dual rank two tensor operator with scaling dimension A = 2 + |m]. 
In particular, this charged operator has angular momentum / = 1 and thus can serve as an 


18T| the theory of superfluidity and superconductivity, a discontinuity of the free energy was discussed 
theoretically and an exactly solvable model for such phase transition was given in ref. [107]. The zeroth 
order transition was also observed in holographic superconductors in refs. [108, 109, 110]. 
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order parameter for p-wave superconductors. Since what we are interested in is a system at 
finite temperature and chemical potential with respect to the global U(1) symmetry, we will 
construct electrically charged asymptotically AdS black holes in gravity side. The normal 
phase with no condensate is described by the electrically charged Reissner-Nordstr6m AdS 
black hole, which is spatially homogeneous and isotropic. This model is specified by two 
parameters m and e. It was shown in ref. [21] that when e? > m?/2 this black hole is 
unstable to developing non-trivial two-form hair that is dual to p-wave superconductors 
with helical order. 


4.3.1 Boundary conditions 


The helical black hole solution was constructed in ref. [22] in which the authors adopted 
the ansatz 


ds? = =g fd? tg dr’ +h w tr’ (et wy + e” w), 


78 
C = (i cı dt + codr) A w2 + C3 w1 A Ws , A=adt, ve) 
where the one-forms w; are given by 
Wy = dx 3 
w = cos (kzı) dzo — sin (kx) dz3, (79) 


w3 = sin (kzı) dz + cos (kzı) dzz. 


Note that the constant t and r slices in the above metric are spatially homogeneous of 
Bianchi type VIlo. All eight functions in the ansatz depend on the radial coordinate r only 
and k is a constant. After substituting the ansatz into the action (76), one finds that cı and 
C can be determined by other functions, thus we are left with six independent functions 
including f, g, h, a, c3 and a. More precisely, f and g satisfy first order differential 
equations and other functions satisfy second order equations. 

To solve the coupled equations of motion for above six functions, one needs to specify 
suitable boundary conditions in the horizon r, and the conformal boundary r — ov. 
Regularity at the horizon demands that g(r;,) = a(r) = 0 and all of them have analytic 
expansion in terms of (r — ra). We then find that the full expansion at the horizon is fixed 
by six parameters, i.e., ra, f (Tn), h(Tn), a(Tn), a’ (Tp) and c3(r;,). Near the boundary r —> oo, 
one demands asymptotically AdS geometry with the fall-off 


g=r? (1—-Mr“+---), f= fol- ear“ +), 
h=r (1+er +), a= car, (80) 
a= fo (u—pr7/2+---), C3 = erT He, 
which is determined by eight parameters M, fo, Ch, Ca; H, P, Cy and k. One should note that 
the expansion of c3 is chosen so that the charged operator dual to the two-form C has no 


source, thus can spontaneously acquire an expectation value proportional to c, which is 
spatially modulated in the x; direction with period 27/k. u and p are regarded as the 
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chemical potential and charge density in the dual system respectively. The holographic 
interpretation of the other UV parameters will be given below. Observe that when k 4 0 
the order parameter rotates in the (£2, x3) plane as one moves along the x, direction thus 
there is a reduced helical symmetry. 

There are two scaling symmetries of the coupled equations which can be used to set 
u = fo = 1. To solve the six differential equations, we need to specify ten integration 
constants. However, we have fourteen parameters in two boundaries minus two for the 
scaling symmetries. Therefore, we expect to leave with a two parameter family of black 
hole solutions which can be selected as temperature T and wave number k. 


4.3.2 Thermodynamics 


We shall work in grand canonical ensemble with the chemical potential u fixed. The 
thermodynamic potential of the boundary thermal state is identified with temperature 
T times the on-shell bulk action in Euclidean signature. We denote w as the density of 
thermodynamic potential per spatial volume in dual field theory. Then one can obtain the 
following expression for the free energy density 1° 


w=—-M=3M+4+8c,—-—pp-—sT, (81) 


where the entropy density s = 4mr?h(rp,) and fo is set to be one. From above equation one 
can immediately obtain the Smarr-type formula 


4M +8c, —pp—sT=0. (82) 
An on-shell variation of the total action for fixed k gives us the first law 
dw = —sdT — pop, (83) 


and hence w = w(T, p). 
The expectation value of the dual stress-energy tensor is given, after setting fo = 1, by 


Ta = 3M = Sey, Tiga, = M + Bey, 
Trot. = M + 8c, cos (2kx1) , 

Trz = M — 8cq cos (2kz1) , 

Tene = — 8 Ca Sin (2kx1). 


(84) 


Obviously the stress-energy tensor is traceless as a consequence of the underlying conformal 
symmetry. We further extract the energy density € = Ta = 3M + 8c, from which we can 
rewrite w = £ — sT — up and the first law takes the more familiar form d¢ = Tds+pdp. The 
average hydrostatic pressure p is defined as minus the average of the trace of the spatial 
components. We get p = M + 8c;,/3, and hence the system satisfies the thermodynamical 
relation £ + p = T's + pp. 


19For more details about this result, please see ref. [116]. 
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4.3.3 Helical p-wave solutions 


We focus on the specific case with m = 1.7 and e = 1.88 7° and set u = fo = 1. Starting 
from the AdS Reissner-Nordstr6m black hole solution, as the temperature is lowered, the 
first instability appears at Te ~ 0.0265 and k = ke ~ 0.550. Below T}, there is a continuum 
of hairy black hole solutions appearing with different values of k. 

Figure 12 summaries the free energy density w as a function of temperature T and wave 
number k. One can see that all hairy solutions have smaller free energy than the normal 
solutions at the same temperature and the transition to the p-wave preferred branch is 
second order. For a given temperature T < JT, there is a one parameter family of solutions 
specified by k, and the most thermodynamically preferred solution is denoted by the red 
line. One can prove that while the general hairy solutions in figure 12 have c, 4 0, the 
solutions on the red line do have vanishing c, [116]. 


Figure 12: The free energy density as a function of T and k for the case of m = 1.7 and 
e = 1.88. The red curve denotes the thermodynamically favored p-wave superconducting 
phase minimizing the free energy with respect to k at fixed T. The black curves correspond 
to curves with constant T. The blue line is for the case of some domain wall solutions. 
Reprinted with kind permission from ref. [22]. 


The helical superconducting order can be fixed by the condensate c, and wave number 
k, which are shown in figure 13 with respect to T for the red line in figure 12. Near Te 
one can find the critical phenomenon c, ~ 1.7 x 10°T?7 (1 — T/T.)'/”, which is the famous 
mean field behaviour. As the temperature is lowered, the red line moves smoothly down 
to sufficiently low temperature at which k = ko ~ 0.256. In particular, the ground state 
at T = 0 is also spatially modulated. 

The T = 0 limit of hairy solutions approach a smooth domain wall solution which 
interpolates between AdS; in the UV and a new IR fixed point with an anisotropic scaling. 


20The main reason for this choice is to obtain real scaling dimensions. For other values of m,e which 
can avoid complex scaling dimensions will give similar results [116]. 
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This fixed point in the IR reads 


gr, f= ir, h= kho, @=a0, 85 

a = aor”, C3 = hrs (20) 
with K, ho, a, ao, Co and z all constants. This fixed point solution is invariant under the 
anisotropic scaling r > Atr, t > A*t, £23 4 AX23 and zı —> xı. All those constants can 
be determined by the equations of motion.?! As a typical example, choosing m = 1.7 and 
e = 1.88, one can obtain 


z œ 1.65, K ~ 0.995, ho œ 0.993, 


(86) 
Qo ~ —0.380, ag œ 0.265, co œ 3.69. 
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Figure 13: The condensate c, and wave number k as a function of T for the thermody- 
namically preferred hairy solutions. The blue dots represent the quantities for the domain 
wall solutions. Used with permission from ref. [22]. 


The domain wall solutions interpolating between the UV fixed point and the IR fixed 
point can be specified by the wave number k [22]. One can see in figure 12 that the T —> 0 
limit of the hairy solutions approach these domain wall solutions (the blue line). Similarly, 
in figure 13, the condensate c, and wave number k = kọ œ 0.256 for the domain wall 
denoted by blue dots smoothly connect with the corresponding black hole solution. 

To summarize, a holographic p-wave model with helical superconducting order is intro- 
duced in this subsection. As the temperature is lowered, a helical superconducting state 
emerges spontaneously breaking both the abelian symmetry and the three-dimensional 
spatial Euclidean symmetry down to Bianchi VII[p symmetry. These homogeneous, but 
anisotropic ground states at T = 0 are holographically described by smooth domain wall 
solutions, which exhibit zero entropy density and an emergent scaling symmetry in the far 
IR. 


21Note that one can set fy = k = 1 by scaling t and z1. 
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Further nature of the model (76) has been well studied in ref. [116]. For example, 
some of the p-wave solutions can exhibit the phenomenon of pitch inversion’? and the 
symmetry of the black hole solutions is enhanced at the pitch inversion temperature. The 
superconducting phase can also be (p+ip) order. Depending on the mass and charge of the 
two-form, both the p-wave and the (p + ip)-wave can be thermodynamically favored. The 
two kinds of orders will compete with each other and there can be first order transition 
between them. 


5 Holographic D-wave Models 


It is remarkable to see that rather simple and generic gravity models can capture many 
features of the phase structure of superconducting systems. Nevertheless, in order to con- 
struct more sophisticated and more realistic models one clearly needs to include additional 
ingredients. The focus of this part is on realising an important missing phase, i.e. d-wave 
superconductivity (superfluidity). The importance is self-evident since many unconven- 
tional superconductors admit either d-wave or mixed symmetry. A natural candidate for 
modelling the d-wave condensate is to use a charged spin two field in the bulk, instead of 
a charged scalar field or a vector field. Based on this approach, there are two acceptable 
holographic models describing the d-wave condensate in the literature. 

The authors of ref. [14] first constructed a minimal gravitational model by introducing 
a symmetric, traceless rank-two tensor field minimally coupled to a U(1) gauge field in 
the background of an AdS black hole. The d-wave condensate appears below a critical 
temperature via a second order phase transition, resulting in an isotropic superconducting 
phase but no hard gap for its optical conductivity. Let us call it CKMWY d-wave model in 
terms of the initials of the five authors. The other effective holographic d-wave model was 
proposed soon after the first one with the same matter fields but with much more complex 
interactions [15]. The phase diagram, optical conductivity, as well as fermionic spectral 
function were investigated in detail. With a fixed metric, this model has advantages such 
as being ghost-free and having the right propagating degrees of freedom. This model will 
be named as BHRY d-wave model for short in what follows. 


5.1 The CKMWY d-wave model 


To construct a holographic d-wave model, the minimal effective bulk action including 
gravity, U(1) gauge field and tensor field reads [14] 


1 6 t v î v 1 v 
S= aa | tv [e+ T2 B (Da Bvy) D" B” — m BB" g gem" ? (87) 


22 As the temperature is lowered, the pitch (27/k) first increases, becoming divergent (i.e., k = 0) at 
some particular temperature, then changes sign and finally decreases in magnitude to a value k < 0 at 
T =0. 
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where D, = V, — iqA, is the covariant derivative in the black hole background, L is the 
AdS radius that will be set to unity, and q and m? are the charge and mass squared of 
Biv, respectively. Working in the probe limit, ie. q — oo with qA, and qB,, fixed, 
the matter part can be treated as perturbations in the 3+1 dimensional AdS black hole 
background (25). 

We would like to realize a d-wave superconductor on the boundary such that a conden- 
sate emerges on the x — y plane with translation invariance and the rotational symmetry 
is broken down to Z(2) with the condensate changing its sign under a 7/2 rotation on the 
x — y plane. Therefore, we use an ansatz for B,, and A, as 


Brr = —Byy _ y(r) , A= olr) dt ’ (88) 


with all other field components being turned off and y(r) and ¢(r) being real functions. 
The background geometry is fixed as AdS-Schwarzschild black hole given in (25). Then 
the final equations of motion read 


1 1 2 2 42 
wd -2w- (H -7 TE) uno, 


a P (89) 
n ! q _ 
gaid- Teo =o. 


These two equations are very similar as the case for the Abelian-Higgs model (see equa- 
tions (27)). Therefore, it is natural to expect ~ to condense spontaneously below a critical 
temperature. More precisely, we demand the following asymptotic form near the AdS 
boundary r — oo 


@=p—p/rt+-:, w= firt, (90) 


with AL = Lawl tamne Note that the expansion of w is chosen such that the charged 
operator dual to B,,, has no source, thus can acquire an expectation value proportional 
to fı spontaneously. According to holographic dictionary, u is interpreted as the chemical 
potential, and p as the charge density in the dual theory. The order parameter of the 
boundary theory can be obtained by reading the asymptotic behaviour of B, i.e. 


ZERI (91 


where (i,j) are the indexes in the boundary coordinates (x,y). In what follows, we shall 
keep the chemical potential u fixed and choose q to be minus one, which is the setup 
adopted by ref. [14]. 

The d-wave condensate as a function of temperature can be obtained numerically, which 
is shown in figure 14. One can see clearly that below T}, the tensor field is Higgsed to break 
the U(1) symmetry spontaneously in the boundary theory. Numerical calculation further 
ensures that the phase transition characterized by the d-wave condensate is second order 
with the mean field critical behaviour fı ~ (T, — T)!/?. The conductivity has also been 
studied, which uncovered that the AC conductivity is isotropic and below T}, the DC 
conductivity becomes infinite but has no hard gap. 
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Figure 14: The d-wave condensate as a function of temperature. The condensate goes to 
zero at the critical temperature Tą. We choose m? = —1/4. Used with permission from 
ref. [14]. 


5.2 The BHRY d-wave model 


The approach of the CKMWY d-wave model just writes down a minimal action for the spin 
two field without looking in detail at the constraint equations required to get the correct 
number of propagating degrees of freedom. Soon, the authors of ref. [15] analyzed in more 
detail the effective action for the spin two field and how the constraint equations could 
be satisfied. The desired theory for a charged, massive spin two field in a fixed Einstein 
background takes the following form 


1 d+1 1 v 
S = sa fa T ~~ gE mE" + La), 
La = -|D ul? + 2D, + Duel? — [Du Dro + h.c] — igh, ~~ (92) 
” 1 
—m? (pul? = ||?) Eg OR ae y à garel ) 


where D, = Vp — iqu, 9 = Pw Po = D* Pup and R*,,, is the Riemann tensor of the 
background metric. The above theory is ghost-free and describes the correct number of 
propagating degrees of freedom. The disadvantage is that one has to be restricted to work 
in a fixed background spacetime that satisfies the Einstein condition R,., = 2 Qu: In 
the context of holographic superconductors, this restriction forces us to work in the probe 
approximation where the spin two field and gauge field do not influence on the metric. One 
such a geometry is given by the AdS-Schwarzschild black hole with a planar horizon 
oe z? d 

ds? = = ( — f(z) dt? + dz? + a fi) =1- (=) (93) 
The black hole horizon is located at z = zp, while the conformal boundary of the spacetime 
is located at z = 0. The temperature of this black hole is 


d 
T= (94) 
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5.2.1 The d-wave condensate 


We consider an ansatz where ,, and A, depend only on the radial coordinate z and only 
the space components of y,,, are turned on. According to ref. [15], it is consistent to turn 
on a single component of Yy and to set other components of the gauge field except for A; 
to be zero. Then our ansatz is 


L? 


A= (2) dt , Pay(Z) a 272 


p(z), (95) 


with all other components of p, set to zero, and ¢ and 4% real. 
With the above ansatz (95), the equations of motion for ¢ and w are given by 


1 d= 2 42 272 
w+ (E-S*) ws (E -T )yno, 


f P f 
3—d eL? (96) 
g+ g-o 
z’ f 


Here the prime denotes the derivative with respect to the radial coordinate z. To solve 
the above coupled equations, one demands that two fields near the boundary z = 0 should 
behave as 


Papp? ter, pp (97) 


where A = d/2 + Vd? + 4m?L?/2. The unitary bound implies that A > d for spin two 
operators. Therefore, the mass of y,,, has a lower bound, i.e. m° > 0.7? Note that the 
fall-off of ~ is chosen so that the dual charged operators have no deformation but can 
acquire expectation value spontaneously. Up to a normalization, the coefficients u, p and 
w+ are interpreted as chemical potential, charge density and the expectation value of the 
xy component for the spin two operator Osy, respectively. At the horizon, one should 
require (z) = 0 in order to keep g” A,A, being finite at the horizon. 

In what follows we will focus on d = 3. The resulting boundary value problem can 
be solved directly, for example, by shooting method. A typical dependence of (Osy) on 
the temperature is shown in figure 15. As we lower the temperature, the normal phase 
becomes unstable to developing tensor hair at a certain critical temperature Te. This is a 
typical second order phase transition. 


5.2.2 Conductivity 


We are interested in the electromagnetic response of the condensed phase. To incorporate 
this feature, we extract the optical conductivity of this d-wave model by linear response 
theory. The conductivity tensor o;; can be defined through 


J; = 017 Ej (98) 


23In fact, the previous d-wave model does not consider this aspect. 


41 


1/A 


<Oxy>""U/p 


2.5 


l 
Ww 
in PALE 


2.0} = 


1.5% 


1.0} 


0.5} 


0.2 04 06 08 ro Te 


Figure 15: The condensate (O,,) in d = 3 dimensional spacetime as a function of the 
temperature for various values of A. The curves form top to down correspond to A = 6,4 
and 3.5, respectively. Used with kind permission from ref. [15]. 


where i,j = x,y. J and E are the electric current and electric field, respectively. To 
compute the conductivity in a holographic framework we turn on a source for the current 
J; dual to the gauge field in the bulk. Following the standard approach discussed in 
previous sections, we perturb the gauge field by 6A = eA, (r)dx. The bulk equations of 
motion couple linear fluctuations of the gauge field A, to some spin two components. To 
obtain a consistent set of equations, one should also turn on time dependent fluctuations 
of the complexified fields Yty, Pla Pzy and ply. The linearized equations of motion for the 


e™™t component of these fluctuations are [15] 


1 2 
# £ W 
jade A Ar e-e etue] 


pe pe oP 
= hy = Fa) + E WF- OF, Pa) > (99a) 
0= y+ T oy — ee Oye? = - we YArs + : [2(w + gd) yi, + 4P Pay] » 
(99b) 
0 = [w+ 40)?2? — mL? F] pay + 7 LPafwAl + ASYA, 
- iw +49) — 5 [Aw + 06)z +082] Gy - (99¢) 


The equations for gh, and ly can be obtained by complex conjugation and an additional 
transformation w to —w from the last two equations. The functions y,, and Ply can be 
eliminated from the first two equations by virtue of (99c), leaving three coupled differen- 
tial equations for Az, Pty and ghy Since the fluctuation A, decouples from above set of 
fluctuations, we can conclude that the Hall conductivity o,,(w) is vanishing. 
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The conductivity is related to the retarded Green’s function for the charge current. 
To calculate the retarded function, one should impose causal boundary conditions on the 
equations of motion. As a consequence, the near-horizon modes of the gauge field and spin 
two field are falling into the horizon, i.e., Az, Yty and yl, have the behaviour as 

(op =a PE, (100) 

Near the boundary z = 0, the asymptotical behaviour for the perturbation fields Az, Pty 
and yh, is given by 

Py = Piy-2^- T Piy+2^t Taa (101) 


A— Ay 
Ply = plz + plyt E 


with A, = HE9 Him It VE kimila Here yy- and T are identified as the source terms, while Yyy,+ 


and yl, , are the normalizable fluctuations. Since the presence of source terms Y- and 
he will also source the U(1) current, one should look for solutions where the source in 


the series expansion of Yy and gl, should vanish. Finally, one can obtain the conductivity 
in the x direction as 


AY 
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Figure 16: The real part of the conductivity as a function of frequency for a A = 4 
condensate. Used with permission from ref. [15]. 


In order to obtain the conductivity in the y direction, we look at the effect of a 7/2 
rotation of the condensate in equations (99). This rotation operation flips the sign of 
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w, which is equivalent to flipping the sign of A, by viewing of the equations of motion. 
Changing the sign of A, is equivalent to flipping the sign of both the electric field (the A® 
term) and the current (the AW term), thus conductivity will be unaffected under such a 
change of sign. This implies that the conductivity is proportional to the identity matrix. 
In fact, the isotropy of conductivity is a consequence of the symmetries that o;; has in the 
d-wave case. An isotropic conductivity for a d-wave superconductor can be also produced 
in an explicit microscopic model [117]. 

The numerical results for the conductivity are shown in figure 16 for a d-wave conden- 
sate of conformal dimension A = 4. As the temperature is lowered one can observe a spike 
in the conductivity, which is a signal of a bound state. This spike is localized at smaller 
values of w as the temperature is lowered. A second spike in the conductivity appears 
to disappear as the temperature is decreased. One can also see that Re[o(w)| does not 
vanish even for arbitrary small frequency w, so there is no hard gap in the dual boundary 
superconducting phase. 

The fermionic spectral function in this holographic superconductor with a d-wave con- 
densate has been well studied in ref. [118]. It was showed that, with a suitable bulk 
Majorana coupling, the Fermi surface is anisotropically gapped. At low temperatures the 
gap shrinks to four nodal points, while at high temperatures the Fermi surface is partially 
gapped generating four Fermi arcs. The (d+id) condensate for the BHRY model was inves- 
tigated in ref. [119], in which the existence of fermi arcs is confirmed and a non-vanishing 
Hall conductivity is obtained in the absence of a magnetic field. 

Although both d-wave models we reviewed above can be used to study the properties 
of a superconducting phase transition with a d-wave order parameter in a dual strongly 
interacting field theory, the construction is not ideal. For example, the BHRY d-wave 
model can only work in the probe limit. However, it is well known that including the back 
reaction of matter fields would lead to a much richer phase structure. It will be desirable to 
study a consistent holographic d-wave model with back reaction. To write down an action 
for a charged spin two field propagating in a curved spacetime is challenging, because it 
usually suffers from non-hyperbolic and non-causal behaviour of the spin two field. Apart 
from those two effective models, the authors of ref. [120] discussed top-down models for 
holographic d-wave superconductors in which the order parameter is a charged spin two 
field in the bulk. 


6 Competition and Coexistence of Superconducting 
Order Parameters 


The holographic models of s-wave, p-wave and d-wave superconductors, which have scalar, 
vector and spin-2 order parameters respectively, have been discussed in the previous sec- 
tions. These models were based on a specific setup where the dynamics in the bulk involves 
only a single order parameter. It is desirable to generalize the single order parameter case 
to multi order parameters case because the real high Te superconducting systems indeed 
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involve various orders, such as magnetic ordering and superconductivity, see figure 17 and, 
e.g., refs. [122, 123]. The holographic correspondence provides us a convenient way to 
investigate the interaction for these orders by simply introducing dual fields in the bulk 
as well as appropriate couplings among them. Following this strategy, several works on 
the competition of multi order parameters in the holographic superconductor models have 
already been made [124, 125, 126, 127, 128, 129, 130, 131, 132, 133, 134]. In the following, 
we will review the competition between two s-wave orders [124], the competition between 
s-wave order and p-wave order [125, 126] and the competition between s-wave order and 
d-wave order [127] one by one. The first case concerns the competition between two orders 
with the same symmetry and the last two cases are to study the competition of orders with 
different symmetry. The phase diagrams are also drawn for the corresponding models. 
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Figure 17: Schematic phase diagrams of the cuprates (left) and the pnictide superconduc- 
tor Ba(Fei—-+Cor)Asə (right). In the right plot, the antiferromagnetic phase is labeled by 
AFM, the normal (tetragonal) phase is denoted by Tet, and the superconducting phase by 
SC. Note the similarities of the phase diagrams. Adapted with permission from ref. [121] . 


6.1 Competition and coexistence of two s-wave orders 


Historically, Ginzburg-Landau theory has proved to be an extraordinarily valuable phe- 
nomenological tool in understanding single-component superconductors. Its generalization 
to the two-component Ginzburg-Landau model (TCGL) was constructed, and its appli- 
cability to the two-band systems studied in refs. [135, 136, 137]. Upon switching on 
the interband coupling between the two components, this model can describe the phe- 
nomenon of the two gaps in materials such as MgBo (s+) [138, 139] and iron pnictides 
(s,_) [140, 141, 142]. Applying the multi-band Ginzburg-Landau theories to the gravity 
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side, the holographic multi-band superconductor model can be realized involving some 
competing scalar fields coupled to a single gauge field. Such system exhibits rich phase 
structure. Next, we will describe the holographic model concretely and show results ex- 
plicitly. 


6.1.1 The holographic model 


Let us start a holographic superconductor model with N scalar hairs in (3+ 1) dimensional 
anti-de Sitter spacetime. The action reads [124] 


N 
S= g | dayan E- FFP +S (-[Vd.—ieg Ade|?—melvel?)—Vaniaals (103) 
k=1 

where ex and my (k = 1, 2, ..., N) are the charge and mass of the scalar field Yg, respectively. 
The term Vinta denotes the possible interaction among bulk matter fields. Here one can 
perform a rescaling of the type A, > 5 Aus Ve > ma to set the charge of the scalar field 
w2 to unity. We are interested in the dynamics and mutual interaction among different 
orders. Here we limit ourselves to the case with N = 2. The concrete model we will study 
is described by the following action 


1 6 1 
S = aa | PvR a = Lm], 
2K? L? e (104) 
1 v 
Lm= — ew BO — [Diy]? — midi? — | Dap? — m3’, 


where we have defined Dip = Vp — i Ay, and D2, = V,,—tA,. The parameter ez controls 
the strength of the back reaction and e;/e is the effective charge of pı or the ratio of two 
scalar charges. 

The hairy black hole solution is assumed to take the following metric form 


dr? 
ds? = — f(r) e O dt? + i 


in) +r“(dxz* + dy’), (105) 


together with homogeneous matter fields 


=r), Y2 =%(r), A= ¢(r)dt. (106) 


The horizon rp, is determined by f(r;,) = 0 and the temperature of the black hole is given 

by 

f'(rp)e x0) 
Ar 

One can use the U(1) gauge symmetry to set pı to be real. After using the r component 

of Maxwell’s equations we can also safely choose pə to be real. We will work in the unites 


T= (107) 
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where L = 1. Then, the independent equations of motion in terms of the above ansatz are 
deduced as follows 


I ' 2 42x 2 
wi q -X + wid as -7 =0, 
1 1 2 2X 2 
usd q —% + Si + (Fe — un = 0, 
1 2 2 2 
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where a prime denotes the derivative with respect to r. 
The gravity background describing the normal phase is just the AdS Reissner-Nordstr6m 
black hole with a planar horizon 
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where rp is the black hole horizon and wu is the chemical potential of the black hole. 


6.1.2 Phase transition 


The two band model is controlled by four model parameters, i.e., m?, m3, e2, and e1/e2. 
Here we will choose m? = 0 and m3 = —2. One may expect that the model admits three 
different superconducting phases. The first superconducting phase corresponds to the case 
with w, # 0 and Y2 = 0 (Phase-I). The second superconducting phase corresponds to the 
case with Y> # 0 and %ı = 0 (Phase-II). The third superconducting phase admits the 
region where both scalars condense simultaneously. 

The numerical results are shown in figure 18. The model admits three kinds of the 
coexisting phase. The first kind is shown in the left panel of figure 18. As we lower 
temperature, the scalar y% first condenses at T, where the superconducting phase transition 
happens; when we continue lowering temperature to a certain value, say Th, the scalar Y2 
begins to condense, while the condensate of yı decreases, resulting in the phase with both 
orders; if one further lowers temperature, the first condensate quickly goes to zero at a 
temperature 73; when temperature is lower than 73, there exists only the condensate of 
pa. This superconducting phase is denoted by Phase-C and it is the case uncovered in the 
probe limit in ref. [129]. The second kind of the coexisting phase is presented in the middle 
plot. It is different from the first one in that the coexisting phase survives even down to 


47 


i 04 
0.4 J 0.4 
0.3 ; 0.3 0.3 


Oy" 0.2 2 (0)"^ 0.2 J (oy! 0.2 


0.1 i 0.1 J 0.1 


0.0 = 0.0 a 0.0 
0.030 0.035 0.040 0.045 0.050 0.055 0.00 0.01 0.02 0.03 0.04 0.010 0.015 0.020 0.025 0.030 0.035 
T T T 


Figure 18: The condensate as a function of temperature for three kinds of the coexisting 
phase. The red curve is for %1, while the blue one is for %2. The condensate for Phase-C 
with e1/e2 = 1.95 and e> = 4 is shown in the left plot. The condensate for Phase-A with 
e,/e2 = 1.95 and e2 = 2 in the middle plot and for Phase-B with e1/e = 1.9 and eg = 1.5 
in the right plot. Three plots were taken from ref. [124]. 


a low temperature. We denote this case by phase-A. Depending on the back reaction, the 
inverse of phase-A is also true: the condensate of pı emerges following the condensate of 
we, and then both orders are always present. This case is labeled as Phase-B drawn in the 
right plot. 

To determine whether those above three coexistence phases are thermodynamically 
favored in their own parameter spaces, one should calculate the free energy of the system 
for each phase. Working in grand canonical ensemble, the chemical potential is fixed. In 
gauge/gravity duality the grand potential Q of the boundary thermal state is identified with 
temperature times the on-shell bulk action with Euclidean signature. The Euclidean action 
must include the Gibbons-Hawking boundary term for a well-defined Dirichlet variational 
principle and further a surface counter term for removing divergence. Note that we are 
considering a stationary problem, the Euclidean action is related to the Minkowski one by 
a minus sign as 


6 1 4 
2K S Euclidean = pes V —g[R + T2 + zz Eml] + | Ër V —h(2K — 7) 
1 Yue a le 7 
tof Pav HCH y + Sg) 
e2 Jr—oo L L 


where h is the determinant of the induced metric on the boundary r — oo, and K is 
the trace of the extrinsic curvature. By using of the equations of motion (108) and the 
asymptotical expansion of matter and metric functions near the AdS boundary, the grand 
potential Q can be expressed as 


2K?NQ 
=e 111 
= (111) 
where V2 = f dxrdy and the constant £ is from the asymptotical expansion of f = r° + 
e/r +--+. For the normal phase given in (109), one has € = —r? — a 


The free energy corresponding to phase-A(C) is drawn in figure 19. From each plot, 
phase-A(C) does have the lowest free energy, indicating that once phase-A(C) appears, it 
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is thermodynamically favored. But for phase-C, there is only a narrow window admitting 
the two orders to coexist. Comparing phase-B with phase-A, the only difference is that Y2 
condenses before Yı. From figure 19, one can also see there exist two transition points in 
both cases. One is the critical superconducting phase transition and the other within the 
superconducting phase, which indicates the fact that our system is multi-band. In order 
to see this clearly, one can define the total condensate as (O) = (O1) A1 + (O24) 4, 
and draw (Q) as a function of temperature in figure 20. As one lowers temperature, (OQ) 
emerges at the critical superconducting phase transition point, then at a certain temper- 
ature inside the superconducting phase it has a sudden increase, where the condensate of 
the other w appears. Such a behaviour is very reminiscent of the one in the real multi-band 
superconductor. 
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Figure 19: The free energy as a function of temperature for Phase-A(C) labeled as the 
solid red curve. The equations of motion also admit three other types of solutions, i.e., 
the normal phase (dotted black curve), Phase-I (dashed green curve) and Phase-II (solid 
blue curve). The curve in the insert of the left plot is the difference of free energy between 
Phase-A and Phase-II. One can see phase-A(C) indeed has the lowest free energy. The 
plots were taken from ref. [124]. 


6.1.3 Conductivity 


In order to ensure the system is indeed in a superconducting state, one should calculate 
the conductivity o. Since now the back reaction is included, one has to consider the 
fluctuations of A, and g. Assuming both perturbations have a time dependence of the 
form e~“*, the final equation of motion to calculate the conductivity can be obtained as 


1 1 2 12 9 2 
Alt E-Da + G— Se Gal + WA = 0 (112) 


Since the conductivity is related to the retarded two-point function of the U(1) current, 
i.e, o = EGP (w, k = 0), one imposes the ingoing boundary condition near the horizon 


Ag = (r — ra) ET [ao + a(r — ra) + a(r — ra)? +-+], (113) 
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Figure 20: The total condensate as a function of temperature for Phase-A. We set e1 /e2 = 
2 and e> = 2. The two special points at T, ~ 0.0488 and T œ 0.0298 correspond to 
the superconducting critical point and the point at which Y begins to emerge in Phase-A, 
respectively. The figure was taken from ref. [124]. 


with ao, a1, az being constants. The gauge field A, near the boundary r — oo falls off as 


AQ) 
A, = AM + — H. (114) 
F 


According to the AdS/CFT dictionary, the retarded Green function can be read as GË = 
1 AM 


ZAO from which one can obtain the conductivity 
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The optical conductivity as a function of frequency in the region with two order parameters 
is presented in figure 21. One can see clearly that the optical conductivity in two band 
model behaves qualitatively similar to the model with only one scalar order discussed in 
figure 2. In addition, from the Kramers-Kronig relations, one can conclude that the real 
part of the conductivity has a Dirac delta function at w = 0 since the imaginary part has 


a pole, i.e., Im[a(w)] ~ L, 


6.1.4 Phase diagram 


Constructing the parameter space is helpful to learn in which region the superconducting 
orders can coexist. One can complete this task by just turning the problem as an eigenvalue 
problem. Focus on the concrete model discussed in this paper, i.e., m? = 0 and m3 = —2, 
a good starting point is to find the critical valve of the ratio e1/e2 such that T is a critical 
temperature at which Yı begins to vanish or emerge. At such a temperature, pı is very 
small and can be treated as a perturbation on the background where only w2 condenses 
1 1 2 m2 e2 2 ex 
Wi -(E-¥ + i + = 2 
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Figure 21: The optical conductivity as a function of frequency at temperature T = 0.0273 
for Phase-A. The red solid line is the real part of the conductivity, while the blue dashed 
line is the imaginary part of the conductivity. Here the parameter e;/e2 = 1.95 and eg = 2 
are taken. There is a delta function at the origin for the real part of the conductivity. 
Figure taken from ref. [124]. 


where {¢, f,x} are functions describing the hairy AdS black hole with only 2 non- 
vanishing. Imposing the appropriate boundary conditions, this equation can be considered 
as an eigenvalue problem with positive eigenvalue e7/e3. The full phase diagram for the 
five superconducting phases is shown in figure 22. 

From figure 22, one has as many as five superconducting phases in the model apart 
from the normal phase. Depending on the model parameters e;/e2 and ez, each phase can 
be most thermodynamically stable in some region of parameter space. As one increases 
the strength of the back reaction, the region for Phase-C with the coexisting behaviour of 
two order parameters only in a narrow window is gradually forced to shrink and finally 
vanishes at e$"’“!, while the regions for Phase-A and Phase-B where both order parameters 
always present enlarge. In this sense, one can conclude that the gravity which provides an 
equivalent attractive interaction between the holographic order parameters tends to make 
the coexistence of two orders much more easy rather than more difficult. 

In this subsection 6.1, a holographic superconductor model with more than one order 
parameter in four dimensions has been studied, where each complex scalar field in the bulk 
is minimally coupled to a same U(1) gauge field. This can be interpreted as a holographic 
multi-band superconductor model. Concretely, we have discussed the two-band case with 
mass squares m? = 0 and m3 = —2 for two bulk scalar fields yı and we, respectively. 
Depending on the strength of the back reaction 1/e3 and the relative charge ratio e1/e2 of 
the two scalar fields, the model admits as many as five different superconducting phases. 
Three of them, denoted by Phase-A, Phase-B and Phase-C, exhibit the coexistence region 
of two order parameters. More specifically, for Phase-C, as one lowers the temperature, 
the second scalar Y% condenses following Yı will completely suppress the condensate of 
the first order, i.e., yı will go to zero finally. The condensate behaviours in Phase-A and 
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Figure 22: The full phase diagram for the five superconducting phases. Depending on 
e,/e€2 and ez, the phase diagram is divided into five parts. The most thermodynamically 
favored phase in each part is labeled. This figure was taken from ref. [124]. 


Phase-B are similar. One of the two orders condenses first, and once the other begins to 
condense, both always coexist. However, this model is limited to the competition of the 
order parameters with the same asymmetry. Therefore it is quite interesting to study the 
holographic models with superconducting order parameters with different spins. This will 
be done in the following subsections. 


6.2 Competition between s-wave and p-wave orders 


In this subsection, we will study two holographic superconductor models with both s-wave 
and p-wave condensed. One is proposed in ref. [125], where the authors built a holographic 
superconductor model with a scalar triplet charged under an SU(2) gauge field in the bulk. 
The other holographic s+p model in ref. [126] consists of a scalar doublet charged under 
an U(2) gauge field living in a planar Schwarzschild black hole geometry. The discussions 
for both models are limited to the probe limit case. 


6.2.1 The holographic s+p superconductor with a scalar triplet charged under 
an SU(2) gauge field 


To realize the s-wave and p-wave superconductivity in one model, we first consider a real 
scalar triplet charged in an SU(2) gauge field in the gravity side. The full action is [125] 


= 1 
OR? 


S f avan — 2A) + Sm, 
(117) 


1 1 
Sm =; f dey g( -DV DIV" — LF F — muwe), 
ge 
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where W* is an SU(2) charged scalar triplet in the vector representation of the SU(2) gauge 
group, and 
D,W* = 8,0? + PAL. (118) 


Fi, is the gauge field strength which is the same as (43) and reads 
Fi, = ð Ag — 0,A% + e” A? AS. (119) 


ge is the Yang-Mills coupling constant as well as the SU(2) charge of Y°. One can redefine 
the fields A’, and W* to get the standard expression where the coupling ge appears in the 
derivative operator D,,. Here we limit ourselves to the case of probe limit. This limit can 
be realized consistently by taking the limit ge —> oo. 

In the probe limit, we consider the d+1 dimensional AdS black brane as the background 
with metric i 
f(r) 


x's are the coordinates of a d — 1 dimensional Euclidean space. The function f(r) is 


par (1-8), a21) 


with r, the horizon radius. Here the AdS radius L has been set to be unity. The temper- 
ature of the black brane is related to r, as 


ds? = —f(r)dt? + dr? + r°dzidz’. (120) 


Pan (122) 
4T 
This is just the temperature of dual field theory in the AdS boundary. 


Let us consider the following ansatz for the matter fields 
Y? = P(r), Al =¢(r), A = Y,(r), (123) 


with all other field components being turned off. In this ansatz, we take A; as the electro- 
magnetic U(1) field. With this ansatz, the equations of motion of matter fields in the AdS 
black brane background read 


wr 4 = | DEZ | os = 0, (124) 
a a a = 0 


One can see Y3? and W, are not coupled in their equations of motion, but they are both 
coupled to the same U(1) electromagnetic field. In this model, thus one can easily realize 
the s-wave and p-wave superconductivity consistently. 
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Figure 23: Condensate of the operators in the s+p coexisting phase. The blue curve is 
for the condensate of the p-wave operator, while the red curve is for the condensate of the 
s-wave operator. The figure was taken from ref. [125]. 


We take the case with A = Aeg = (6+ V3)/4 as an example. The condensate behaviour 
for the coexisting phase is drawn in the left plot of figure 23. We can see that the s + p 
coexisting phase starts from the p-wave phase and ends with the pure s-wave condensate 
phase. Based on the calculation of free energy shown in figure 24, we confirm that the s+p 
coexisting phase indeed has the lowest free energy and is thus thermodynamically favored 
in the temperature region. Thus the potential first order phase transition from the pure 
p-wave phase to the pure s-wave phase is replaced by the phase transitions from the p- 
wave phase to the s-wave phase through an s+p coexisting phase. And all the three phase 
transitions are continuous ones, and are of characteristic of second order phase transition 
within the numerical accuracy. 

From figure 24, we see that the Gibbs free energy curves of the s-wave and p-wave phases 
have an intersection when Aer < A < Aaz. The s+p coexisting phase just exists in this 
interval. By computing the values of T3?' and T? and getting the relations T3?'(A) and 
T? (A) in the region Aer < A < Aur , a phase diagram of the holographic model on the 
A-T plane can be shown in figure 25. We can see from the figure that the system contains 
four kinds of phases known as the normal phase, the s-wave phase, the p-wave phase and 
the s+p coexisting phase. The s+p coexisting phase is favored in the area between the 
blue line and the red line. The region for the s+p coexisting phase is very narrow in the 
phase diagram. This is similar to the situation of the coexisting phase with two s-wave 
orders in the probe limit studied in ref. [129]. However, for the latter, the region with the 
coexisting phase is enlarged with the full back reaction [124]. This would be due to the 
additional interaction between the two scalar fields in the bulk through gravity and this 
interaction reduces the repellency between the two condensates. 

Very recently, the back reaction effect was included in this model [143], which showed a 
rich phase structure and various condensate behaviours such as the “n-type” and “u-type” 
ones. The phase transitions to the p-wave phase or s+p coexisting phase become first 
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Figure 24: (Left) The Gibbs free energy versus temperature for various phases. The 
black solid curve is for the normal phase, the blue solid curve is for the p-wave phase, and 
the dashed lines from bottom to top are for the s-wave phase with operator dimension 
A = Az, (6+ V3)/4, and Aaz, respectively. (Right) The Gibbs free energy in the region 
near the intersection point of the p-wave curve and the s-wave curve with A = (6+ V3) /4. 
The plots were taken from ref. [125]. 


p-wave phase 


Figure 25: The A — T phase diagram. The normal phase, the s-wave phase, d-wave phase 
and the coexisting phase are colored differently. The right plot is an enlarged version of 
the coexisting region. The figures were taken from ref. [125]. 


order in strongly back reacted cases. The phase diagrams similar as figure 25 in different 
strength of back reaction were constructed, indicting that the region for the s+p coexisting 
phase is enlarged with a small or medium back reaction parameter, while is reduced in the 
strongly back reacted case. 


6.2.2 The holographic s+p superconductor with a scalar doublet charged un- 
der a U(2) gauge field 

In this sector, we consider a holographic s+p model consisting of a scalar doublet charged 

under a U(2) gauge field living in a (3 + 1)-dimensional Schwarzschild-AdS black brane 
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geometry. The action for the matter sector reads [126] 


S= | dxy=9 (-37 Fi, -mwt —(DNNID,W) , (125) 


with 
v=v2(%) | Dy=O,-iA,, Ay= ACT:, (126) 
1 1 
To = 9 , = act? 


The system lives in the Schwarzschild-AdS background (120). Considering the following 
consistent ansatz for the fields 


AY =8(r), A? =O), AP =ulr), v=H(r), (127) 


1 _Q\2 2 2 
w+ (E4 jy+ (E or eno, 


f 4f? f 4f 
2 2 
Pp y= e0) = 0, 
E (128) 
ie a =) ee = 
f ry 
f' o2 yp? 
w" + =w + >w- w= 
f f? f 
In what follows we choose the scalar to have m? = —2 and thus the corresponding dual 


operator has mass dimension 2. 
The UV asymptotic behaviour of the fields, corresponding to the solution of equa- 
tions (128) in the limit r — oo, is given by 


=u- pyri 5 
O = p3 — p3/r +, 
TEM OTEN Pon on: 


D fp $y fr? pes 


where in the dual field theory side, u and p are respectively the chemical potential and 
charge density corresponding to the overall U(1) C U(2) generated by To, whereas u3 and p3 
are the chemical potential and charge density corresponding to the U(1) C SU(2) generated 
by T3. py is p aes of a scalar operator of dimension 2, while y® is its expectation 
value. Pinay w and w are the source and vacuum expectation value of the current 
operator J}. Notice that in a background where w(r) condenses the SU(2) c U(2) is 


(129) 
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spontaneously broken, and moreover spatial rotational symmetry is spontaneously broken 
too. 

We are looking for solutions of the equations (128) where Y, w, or both acquire non- 
trivial profiles. First we will switch on a chemical potential u along the overall U(1), while 
requiring that the other chemical potential u3 remains null. Therefore the UV boundary 
conditions are 

yY=0, w=0, ps =0. (130) 


In the IR regularity requires A; to vanish at the black hole horizon. So far, the holographic 
multi-component superfluid model has been realized. 
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Figure 26: Left: Condensates Y® (solid) and w® (dashed) as a function of 1/p in the 
s-wave (blue) and s+p-wave (red) phases. The p condensate appears at Hsp such that 
Ls/ [sp = 0.395. The inset zooms in on the plot of 7) to show the difference in the scalar 
condensate between the s (blue) and the s+p (red) solutions. Right: Free energy of the 
different solutions versus 1/u: normal phase in black, s-wave phase in blue, and s+p-wave 
phase in red. Reprinted with permission from ref. [126]. 


In the left plot of figure 26 the condensates (O2) ~ 7) and (I?) ~ w are plotted 
as a function of the chemical potential. Notice that the solution where both condensates 
coexist extends down to as low 1/y as where the decoupling limit is trustable. And the free 
energy for the different solutions is shown in the right plot of figure 26. At small chemical 
potential only the normal phase solution exists. At u = us ~ 8.127 there is a second order 
phase transition to the s-wave solution. If one keeps increasing ju, at Hsp ~ 20.56 there is a 
second order phase transition from the s-wave phase to the s+p-wave phase. The system 
stays in the s+p-wave phase for u > Hsp. 

Next, we relax the condition u3 = 0 and study the phase diagram of the system as 
a function of u and u3/u. Notice that turning on a second chemical potential means to 
explicitly break U(2) + U(1) x U(1). The system can now be interpreted as a holographic 
dual to an unbalanced mixture. The UV boundary conditions now read 


ge =0,, wav (131) 
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Figure 27: Phase diagram of the unbalanced system as a function of 1/u and u3/u. Second 
order phase transitions are denoted by blue lines, whereas the red line corresponds to a 
first order phase transition. Reprinted with permission from ref. [126]. 


In this case, by computing the free energy of the different solutions, the phase diagram 
of the system as a function of 1/j: and ju3/p is plotted in figure 27. For small values of ju3/11, 
the solution where both condensates coexisting extends down to as lower 1/j as where we 
can trust the decoupling limit. As |3|/j gets larger, the transition to the s+p-wave phase 
happens at a higher value of u. For |ju3|/j is large enough, the p-wave phase is preferred at 
intermediate values of u. Therefore, as u is increased above a critical value up the system 
goes from the normal to the p-wave phase through a second order phase transition. If u is 
increased even further a first order phase transition takes the system from the p-wave to 
the s-wave phase. The tricritical point where the normal, s-wave and p-wave phases meet 
happens at 1/p œ 0.223 and |3|/u ~ 0.815, the p-wave solution is never energetically 
preferred for |us|/u < 0.815. 

In summary, in this subsection we have reviewed the competition between s-wave and 
p-wave order through two holographic superconductor models. The first model is realized 
the s+p superconductor with a scalar triplet charged under an SU(2) gauge field and the 
other is constructed with a scalar doublet charged under a U(2) gauge field. The s+p 
coexisting phase exists in both models. In the first model, the s+p coexisting phase is 
narrow and one condensation tends to kill the other. This competing behaviour is similar 
to the case shown in the condensed matter system [144]. However, in the second case, the 
condensates feed on different charge densities and the coexisting phase survives down to 
a low temperature. Therefore, it should be noted that the competing scenario is model 
dependent. In next subsection, we will study the competition between the scalar field and 
the tensor field, i.e., the competition between the s-wave and d-wave orders. 
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6.3 Competition between s-wave order and d-wave order 


In section 5, we have mentioned two acceptable holographic models describing the d-wave 
condensation. The CKMWY d-wave model is reviewed in subsection 5.1 and the BHRY 
d-wave model in subsection 5.2. In order to realize the condensation of s-wave order and 
d-wave order in one holographic model, one can simply combine the Abelian-Higgs model 
with a d-wave model. Thus, one could have two holographic models with s-wave order 
and d-wave order. Here we will discuss the competition between s-wave order and d-wave 
order for both d-wave models in the probe limit where one neglects the back reaction of 
matter fields to the background geometry [127]. The phase structures are given and the 
behaviours of the thermodynamic quantities for the s+d coexisting phase are also studied. 
The coexisting phase does appear in both models and is thermodynamically favored. 


6.3.1 The s-wave + BHRY d-wave model 


To study the competition between s-wave and d-wave orders, let us first start with the 
holographic model by combining the Abelian-Higgs s-wave model and BHRY d-wave model. 
The holographic model with a scalar field %1, a symmetric tensor field y,,, and a U(1) gauge 
field A, is described by the following action [127]: *4 


1 1 
S= 55 | daya EE" — Dia -måla + La, 
La = -|De + 2|D,e"? + Duel- [Due Dup + h.c] — ig Fp yy (132) 


P 1 
=mi (pul? = ll’) + Rupee à grel, 


where D, = V, — iqA, and D, = V,—imA, P = Pm Yo = g Dy Pup and FP yas 
is the Riemann tensor of the background metric. %1 is the scalar order and Wy, is the 
tensor order. The parameters qı and qz are the charges of the scalar and the tensor fields, 
respectively. One can perform a rescaling to set the charge qı of the scalar to be unity. 
Then the phase structure of this theory is determined by the ratio q2/q1 by fixing the mass 
square of the scalar field m? and the mass square of the tensor field m3. We shall set qı = 1 
without loss of generality in the following discussion. 

Working in the probe limit, we choose the background metric to be the 3+1 dimensional 
AdS-Schwarzschild black hole with planar horizon (120). And we consider the following 
ansatz 


Apu dx" = olr) dt , V1 = (r) Pry = Pyr = F palr) ’ (133) 


2 


?4Ref. [128] also discussed the following model, but with a coupling between the scalar field and the 
tensor field, and studied the phase structure in terms of the coupling parameter and temperature with 
fixed charges of two orders. In the folloowing discussion, there is no direct interaction between scalar and 
tensor fields and the model parameter is the ratio of two fields. Note that in paper [128], when the coupling 
n = 0, there also exists a coexisting phase under the model parameters m? = —2,m3 = 0 and q2 = 1.95. 
Both results are consistent with each other in that case. 
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with @(r), yı(r) and we(r) all real functions. 
With the above ansatz (133), the equations of motion for ¢@, yı and wy» are given by 


9 1 2 2 
gE- yyt- Sout =o, 
1 2 2 2 
(+ Ea + Ey + Gun = hy, =0, (134) 
1 2 242 2 
wh 4 rus th 4 2 ho » =0 


Here the prime denotes the derivative with respect to r. With this model at hand, we can 
study the competition mechanism between the s-wave order and d-wave order. It is easy 
to see that equations (134) have a symmetry 


mi m3, g > 1/q2, 6 > Go, Vi > G22/V2, P2 > V2. (135) 


Under this symmetry transformation, the role of s-wave and d-wave would interchange 
each other. Without loss of generality, here we focus on the case m? < må. 

Before solving the set of coupled equations (134) numerically, we make a briefly qual- 
itative analysis on the possible phases for such a model. Following ref. [129], we rephrase 
the equations for the s-wave and d-wave as a potential problem. The evolution equations 
for s-wave and d-wave in equations (134) can be rewritten as follows 


2 


1 = Viers(y)vr = 0, 


r (136) 
CO E E 
aye? reff (y)W2 
where dy = =f with z = 1/r and Vies(z) = -P3 — a + £2 z3) and Voerr(z) = 
2 12 2 
= P = r | le 2°). Now in terms of the new variable y, the equations of motion 


for s-wave and d-wave are rephrased as a potential problem on a semi infinite line, i.e., 
y € [0,co). Our qualitative discussion is based upon the lemma proven in ref. [129]. For 
the case q? < 1, no matter which gauge field configuration we choose, we always have 
Vieff < Vaerr- Therefore the phase structure of the system is the same as that of s-wave 
holographic superconductor with a single scalar. While, for the case q2 > 1, one may expect 
that the d-wave field with large charge q2 will always dominate. However, the potential 
Vier diverges like 4 near the boundary y = 0 when we lower the temperature. Therefore, 
lowering the temperature possibly makes the mass dependent potential more important 
and hence the s-wave tends to dominate. We will confirm this with the following numerical 
calculation. 

Here we set the mass square m? = —2 and m2 = 7/4 and we take q2 = 2.66 as a 
typical example. Our numerical results confirm that the model does admit the coexistence 
region of two orders with different symmetry, which is drawn in figure 28. We find that 
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the s+d coexisting phase starts from the d-wave phase and ends with the pure s-wave 
condensate phase. The calculation of the free energy confirms that the coexisting s+d 
phase is thermodynamically favored as shown in figure 29. 


Figure 28: Condensate of the operators in the s+d coexisting phase. The blue curve is 
for the condensate of the d-wave operator, while the red curve is for the s-wave operator. 
The figure was taken from ref. [127]. 
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Figure 29: The left plot shows the difference of Gibbs free energy between the supercon- 
ducting phase and the normal phase. The blue curve is for the d-wave phase, the green 
line is for the s-wave phase, while the red curve is for the s+d coexisting phase. The right 
plot is an enlarged version of the left one to show the s+d phase more clearly. Figures 
taken from ref. [127]. 


Based on the above discussions, it can be seen clearly that there exist three particular 
points at which the derivative of the charge density with respect to temperature is discon- 
tinuous, indicating a second order phase transition. The one with the highest temperature 
is the critical point for the superconducting phase transition, while the remaining two 
points are inside the superconducting phase, indicating the appearance and disappearance 
of coexisting phase. We can also see the signal of phase transition from the behaviour of 
the total charge density as a function of temperature and the ratio p,/p versus temperature 
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shown in figure 30, where ps is the superconducting charge density ps = p — pn and pn is 
the normal charge density carried by the black hole. 
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Figure 30: Left: The total charge density as a function of temperature. Right: The ratio of 
the superconducting charge density over the total charge density p,/p versus temperature. 
The plots were taken from ref. [127]. 


In order to ensure the system is indeed in a superconducting state, we calculate the 
optical conductivity o(w), which corresponds to the red line shown in figure 31. We see 
that much more interesting phenomena happen in the low frequency region. Unlike the 
s-wave case which only has a bump at w/T ~ 400 in figure 31, for pure d-wave condensate, 
apart from a much more obvious bump at w/T œ 500, Re(oz.) has an additional spike at 
a lower frequency. This spike may indicate the existence of a bound state [15]. One can 
see clearly that such peak becomes much more sharper in the s+d coexisting state, thus 
the bound state is enhanced due to the additional condensate of s-wave order. 
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Figure 31: The real part (left) and imaginary part (right) of the conductivity as a function 
of frequency at temperature T = 0.018. The red curve is for the s+d coexisting phase, 
the green line is for the pure s-wave phase and the blue curve for the pure d-wave phase. 
Plots taken from ref. [127]. 


To be complete, we give the phase diagram with m? = —2 and m3 = 7/4 shown in 
figure 32, which can tell us in which region the coexisting phase appears. From figure 32, 
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Figure 32: The q2-T' phase diagram. The four phases are colored differently and we label 
the most thermodynamically favoured phase in each region. The figure was taken from 
ref. [127]. 


we see that the coexisting phase exists only in a narrow region in the phase diagram. We 
denote the critical temperature for a single s-wave or d-wave starting to condense as Tes 
and Tea. If we set the charges of the s-wave and d-wave fields to unity, then Tes/u ~ 0.0588 
and Tea/u œ 0.0253. We see that 


e In the regime q2 < Tes /Tea ~ 2.323, the s-wave dominates the system and there is no 
condensation of the d-wave order. 


e As q increases beyond 2.323, the s+d phase appears, which emerges from the d-wave 
phase and ends with a pure s-wave. 


e If we continue increasing q to the case q > 1.155T.,/Tig ~ 2.683, the s-wave order 
never condenses and the resulting phase diagram is the same as that of model with 
only d-wave order. 


Finally, we try to give a qualitative explanation on the mechanism through which the 
condensation of one order affects the dynamics of the other order. Note that here the back 
reaction is not taken into account. Thus the two fields interact only through their effect 
on the gauge field once one or both has (have) condensed. Through looking at the gauge 
field we may give some insight into the competing mechanics between two orders. 


e First, after the d-wave order condenses, if one keeps lowering the temperature and 
reaches the critical temperature at which the pure s-wave would condense, this con- 
densation does not happen. This is due to the fact that the condensation of the 
d-wave increases the effective mass of the s-wave, thus prevents the instability of the 
s-wave to happen, which can be seen from figure 33. This reflects the competition 
between s-wave and d-wave. 
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Figure 33: The blue curve is the effective mass square of s-wave without the condensation 
of d-wave. The red curve is the effective mass square of the s-wave under the condensation 
of d-wave. It can be seen clearly that the effective mass of s-wave increases after the 
condensation of d-wave. This figure was taken from ref. [127]. 


e However, further decreasing the temperature, the condensation of s-wave does hap- 
pen. This is due to the fact that the effective mass of the s-wave is lowered and 
ultimately even if the condensation of the d-wave depleted the gauge potential, the 
background with only d-wave order becomes unstable. 


e Finally, the condensate of the s-wave order kills the first one. This should be due to 
the effective mass square of the s-wave being lower. 


It should be noted that this phenomenon is model dependent. This narrow coexistence 
region of two superconducting orders and the fact that one condensate can eventually kill 
the other also happen for two s-wave orders in ref. [129] and the s+p case in ref. [125]. 
The competition diagram is similar to the competition between the conventional s-wave 
and the triplet Balian-Werthamer or the B-phase pairings in the doped three dimensional 
narrow gap semiconductors, such as Cu,BigSe3 and Sn,_,In,Te in the condensed matter 
system [144]. Although in ref. [144] the competition is apparently between a s-wave order 
and a p-wave order, d-wave and p-wave are similar in some circumstances, for example, 
their excitations of the normal component can be probed using low frequency photons. 


6.3.2 The s-wave + CKMWY d-wave model 


With the same strategy, in this subsection we study the competition between s-wave order 
and d-wave order in the model combining the Abelian-Higgs s-wave model [11] with the 
CKMWY d-wave model [14]. The full action including a U(1) gauge field A,,, a complex 
scalar field yı and a symmetric, traceless tensor field B,, takes the following form [127] 


1 1 7 
OS ae [ee =9(— FFP — |Dyi|* — mill’ + La), (137) 
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with g - - 
La = -g (D Ba) DaB” — mB}, B”. (138) 
Here D, = Va — iq1A, and D, = V,, — 1q2A,. In the probe limit, matter fields can be 


treated as perturbations in the 3+1 dimensional AdS black hole background (120). Let us 
consider the following ansatz 


pi = (r), Boy = —Byy — plr), A, = o(r)dt, (139) 


with all other field components being turned off and w(r), Yə(r) and ¢(r) being real 
functions. Then the explicit equations of motion read 
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Figure 34: The let plot shows the condensation in the s+d coexisting phase. The right 
plot shows the differences of Gibbs free energy between superconducting phases and the 
normal phase. Here the blue line stands for the d-wave phase, the green one for the s- 
wave phase and the red one for the s+d coexisting phase. The figure is reproduced from 
ref. [127]. 


The numerical results are shown in figure 34 for the case with qu = 1, q = 1.34, 
m? = —2 and m3 = —13/4. As seen, a new phase with both s-wave order and d-wave order 
coexistence can appear near 7°"°** and this s+d coexisting phase has the lowest free energy 
and is thus thermodynamically preferred to the s-wave phase and d-wave phase. In more 
detail, as we lower the temperature of the system, it first undergoes a phase transition from 
the normal phase to the pure d-wave phase at T?. Then at T$% a new phase transition 


occurs, and the system goes into an s+d coexisting phase. Finally the system undergoes 
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the third phase transition from the s+d coexisting phase to a pure s-wave phase at TS”. 
Note that all the three phase transitions are second order. 

The feature of the phase transitions can also be seen clearly from the charge density as 
the function of temperature in figure 35. One can see that the charge density with respect 
to temperature is continuous, but its derivative is discontinuous at three special points, 
indicating three second order phase transitions. These features are the same as those for 
the model in the previous subsection. But there is a little difference in the behaviour 
of the total charge density for the d-wave phase. In the s-wave + BHRY d-wave model, 
the total charge density changes monotonously with the temperature, while it behaves 
non-monotonous in the present case. 
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Figure 35: The total charge density as a function of the temperature. The red curve is for 
the normal phase, while the blue one corresponds to the superconducting phase. There are 
three special temperatures at which the derivatives of charge density with the temperature 
are discontinuous. Figure taken from ref. [127]. 


The information of the phase transitions can also be revealed via the behaviour of the 
ratio p;/p with respect to the temperature. From figure 36, one can see that the ratio 
ps/p also has a small kink in the region of the coexisting phase. Comparing figure 30 
with figure 36, we see that in the former case, the green dashed curve for the pure s-wave 
phase intersects with the blue dashed curve for the pure d-wave phase. In contrast, the 
green dashed curve in figure 36 is always lower than the blue dashed curve. Therefore, as 
one lowers the temperature, the ratio ps/p in the s+d coexisting phase increases for the 
former (132), while it decreases for the latter (137). The authors of ref. [125] investigated 
an s+p coexisting phase and found the decrease of the ratio p,/p in the coexisting phase, 
similar to figure 36. They suggested that it might be an experimental signal of the phase 
transition from a single condensate phase to a coexisting phase. Nevertheless, the results 
here uncover that the ratio p,/p versus temperature is model dependent. 

The phase diagram for the model (137) with m? = —2 and m3 = —13/4 in the q2 — T 
plane is shown in figure 37. As the s-wave + BHRY d-wave model, the system also 
contains four kinds of phases known as the normal phase, s-wave phase, d-wave phase and 
s+d coexisting phase. The normal phase dominates in the high temperature region, the 
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Figure 36: The ratio of the superconducting charge density over the total charge density, 
p;/p, with respect to the temperature. The red curve describes the ratio p,/p when the 
system transfers from the d-wave phase to s-wave phase through the s+d coexisting phase. 
The green dashed blue curve is for the ratio p,/p of the pure s-wave phase and the blue 
dashed curve is the ratio for the pure d-wave phase. The figure was taken from ref. [127]. 


s-wave phase dominates in the lower temperature region with small gz below the red curve, 
and the d-wave phase dominates in the higher temperature zone with large qə above the 
blue curve. The s+d coexisting phase is favoured in the area between the red and blue 
curves. The region for the s+d coexisting phase is very narrow in the phase diagram, 
which indicates that the s-wave and d-wave phases generally repel each other, but they 
can coexist in a very small range of temperature. 
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Figure 37: The q»-T phase diagram with m? = —2 and m3 = —8. The most thermody- 
namically favored phase in each part is labeled. The s+d coexisting phase exists only in a 
narrow region. The right plot is an enlarged version for the coexisting region in order to 
see this more clearly. Plots taken from ref. [127]. 


Comparing the two holographic setups, i.e., the model (132) and the model (137), one 
can see some common features as follows. 
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e The s+d coexisting phase does exist in a region of the model parameter q2/q1. Once 
the coexisting phase appears, it is always thermodynamically favoured, compared to 
the pure s-wave and pure d-wave superconducting phases, which can be seen from 
the free energy in figure 29 and figure 34. 


e All phase transitions are second order in these two holographic models. 


e One can see from figure 32 and figure 37 that the phase structure is very similar for 
both models. The region for the s+d coexisting phase is very narrow in the phase 
diagram, indicating that the s-wave and d-wave phases generally repel each other. 


There exist also some differences in the two models. For suitable model parameters 
in the first model, as the temperature is lowered, the s-wave order condenses inside the 
d-wave order resulting in the coexisting phase with both orders. However, when the scalar 
order condenses the first one starts to disappear, and finally only the s-wave condensate 
is left for sufficiently low temperatures. If one changes the model parameter m? + m3, 
the inverse is also true: the condensate of d-wave order emerges following the condensate 
of s-wave order, and then the d-wave condensate finally kills the s-wave order. Those two 
kinds of coexisting phase are one to one correspondence. In contrast, in the second model, 
one sees only the first kind of the coexisting phase. What’s more, for the first model, the 
ratio p,/p increases in the s+d coexisting phase as the temperature is lowered, while it 
decreases in the second case. This gives an obvious evidence that the ratio p,/p versus 
temperature is model dependent. 


7 Coexistence and Competition of Magnetism and 
Superconductivity 


The novel paired mechanism makes p-wave superconductor have many features which are 
different from the traditional knowledge coming from s-wave superconductor both in the- 
ories and experiments. In the usual picture, superconductivity and magnetism are in- 
compatible with each other. Especially, ferromagnetic phase, a spontaneously magnetized 
phase which has nonzero magnetic moment without external magnetic field and appears 
when the temperature is lower than a critical one called “Curie temperature”, can not 
coexist with superconductivity at a sample.” This is rooted in the microscopic theory 
of superconductivity from BCS theory. However, this understanding is broken by p-wave 
superconductor. The discovery of the superconducting ferromagnet?’ materials, such as 
UGep [145], URhGe [146], UCoGe [147] and ZrZn [148], came as a big surprise. In this 


25 However, under special conditions superconductivity may coexist with antiferromagnetic order, where 
neighboring electron spins arrange in an antiparallel configuration. Since antiferromagnets don’t have net 
magnetism, we won’t involve them here. 

26We will use “superconducting ferromagnet” to denote the materials whose Curie temperature is higher 
than superconducting transition temperature and “ferromagnetic superconductor” to denote the opposite 
case. 
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material, superconductivity is realized well below the Curie temperature, without expelling 
the ferromagnetic order. 

The nature of superconducting state in ferromagnetic materials is currently under de- 
bate. For a review of phenomenological theory of ferromagnetic unconventional super- 
conductors with spin-triplet Cooper pairing of electrons, one can see refs. [149, 150, 151]. 
However, the microscopic theory about the coexistence of magnetism and superconductiv- 
ity in strongly interacting heavy electrons is either too complex or insufficiently developed 
to describe the complicated behaviour. So it is still a fascinating thing to find a suitable 
theory to describe the coexistence and competition of the ferromagnetism and supercon- 
ductivity in strong correlated system. 

Holographic frame to discuss the coexistence and competition between spontaneously 
magnetic order phase and superconductivity initiated in refs. [132, 152]. Because of lack- 
ing an individual model to describe spontaneously magnetization and the time reversal 
symmetry broken, these models cannot give complete features of this topic. A very new 
idea proposed in ref. [153] tries to give an independent model describing spontaneously 
magnetization in holographic frame, which opens a new direction. We will introduce the 
main results in this framework. For more details, one can refer to refs. [153, 154, 155]. 


7.1 The holographic model for ferromagnetism/paramagnetism 
phase transition 


Before going on the topic of coexistence of ferromagnetism and superconductivity, let’s 
first review how to build a holographic ferromagnetism which is independent on supercon- 
ductivity in ref. [153]. This model is realized by adding a real antisymmetric field into 
Einstein-Maxwell theory in a (3+1) dimensional AdS spacetime, 


1 6 v 
ARa | tavar aa ay X Lm) (141) 
where 
ee ie eer J 
Lm = = M Vp Mor — a My ~~ a Piy = gV (Mw). (142) 


Here 2k? = 167G and G is the Newtonian gravitational constant, \ and J are two constants, 
m is the mass of the real tensor field M,,,, A, is the gauge potential of U(1) gauge field. 
The antisymmetric tensor field M,,, is the effective polarization tensor of the U(1) gauge 
field strength F,,, with the self-interaction V(M,,,,) which should be expanded as the even 
power of M,,,. The probe limit corresponds to A — 0. Under this limit, the equation for 
polarization field decouples from the gauge field and gravity field, 


VM — mM, — JM, M5 M,, — Fe =0, (143) 
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with the dyonic Reissner-Nordstr6m (RN) background [156], 
r?(—f (r)dt? + dr? + dy’) + a 
rf(ry’ 


| (144) 


ds” 
f(r) = 1- r3 +Ë r4 
A, = p(l—1/r)dt+ Brdy. 


Here the horizon radius has been scaled to be unitary. If we only care about the magnetic 
part of polarization field, then a self-consistent ansatz for polarization field is M,, = 
—p(r)dt \ dr + p(r)dx ^ dy. Taking this ansatz into equation (143), we have 
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where a prime denotes the derivative with respect to r. It is interesting to see that these 
two equations decouple from each other in this case, which makes it to be possible that 
we can neglect the dynamic of p(r) if we only care about the dynamic of magnetism. At 
the horizon, we need to impose a regular boundary condition. Near AdS boundary, the 
linearized equations have following asymptotic solutions for p(r), 


p” + 


(145) 


P ~ par J4 m? 


with 6 = V17 + 4m?. In order to make the theory self-consistent and spontaneous conden- 
sation appear, we need following restriction on parameters, 


3 
-4< m < =y and p} =Q. (147) 


According to the action (141) one can derive the magnetic moment from polarization field, 


which reads PE 

N=-5f dre. (148) 

2 Ji r? 

Here we have set the constant A = 1 in this expression for convenience. This integration 
converges only when p, = 0. In the case without external magnetic field, i.e. B = 0, if 
there is a solution such that p(r) 4 0, the magnetic moment then is nonzero, which gives 
a ferromagnetic phase for dual boundary. Because the action (141) implies transformation 
for p(r) such as p(r) + —p(r) under the time reversal transformation, the condensed phase 
of p gives a time reversal symmetry broken spontaneously, which is necessary for magnetic 
ordered phase. 
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7.2  Ferromagnetism and p-wave superconductivity 


Once the two independent models for ferromagnetic phase transition and p-wave super- 
conductor are in hand, we can combine them to discuss the possibility of coexistence. For 
example, we can combine the Einstein-Maxwell-complex vector theory for p-wave super- 
conductor with ferromagnetic model. The complete action reads [154] 


6 
S= fiv i oes X (L, + Lu + Lom), (149) 
with 
1 V £ V 
Ly = —5Plw ol” — mipko" + iqypup F — Vp, 
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Here a # 0 and O are two coupling constant. Lwm is the Lagrangian for polarization field 
which is just as the same as (142). L, is the Lagrangian for complex vector field, which 
is similar to the one we discussed before. However, there is an additional term V, which 
describes the magnetic moment interaction of complex vector field. This term is irrelevant 
for the previous section where we only care about superconductivity but is relevant when 
we care about spontaneous magnetization. 

Under the probe limit À > 0, a self-consistent ansatz of action (149) is, 


Mw = —p(r)dt Adr + h(r)dx Ady, py = prdx + ipydy. (151) 


Then we can get the equations of motion for complex vector field and polarization field 
under the background (144), 
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where we have defined c(r) as p,(r) = c(r)p,(r). Note that the equation for p(r) decouples 
from the others. The linearized equations near the AdS boundary give following asymptotic 
solutions 7 

Px = m + pa r T2, c= cyr” FE 


h(r) = hart)? 4 hr =8)/2, (153) 


27The asymptotic solution of c(r) depends on the source free condition of pz. When pz+ 4 0, asymptotic 
solution of c(r) becomes c = c} + c_r~*. 
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where 6; = \/1+4m? and 62 = y/17 + 4m2 with m? > —1/4, m3 > —4. As the pre- 
vious subsection, we should impose the condition h} = 0 for the polarization field and 
Px = Cz = 0 for the complex vector field, i.e., we require that the condensation and 
magnetization would happen spontaneously. The equations have solutions only when 
c(r) = 0,41. Because of the equivalent of a —> —a and c + —c, we assume a > 0 
without loss of generality. The magnetic moment is defined as the same as (148), 


oh 
N= -f dr. (154) 


According to the dictionary of AdS/CFT, the expectation value of p-wave superconducting 
order parameter is a complex vector P, whose mode is P = v1 + c?|ps-|. Though the 
expression of magnetic moment density does’t contain the terms of complex vector field, 
it is effected by p,, through the mixture terms in equations (152). 

In the pure p-wave model, the global U(1) and spatial rotation symmetries are broken 
spontaneously when py or py is nonzero without source. Here it is also true. Moreover, 
there is an additional symmetry breaking. If one notes following rules for time reversal 
transformation, 

hah. Py > Ny (155) 


then when h ¥ 0 or py = +p, # 0 (they both lead nonzero magnetic moment), the time 
reversal symmetry is broken spontaneously, which agrees with the fact that a spontaneously 
magnetized phase is with a time reversal symmetry broken spontaneously. 

Because the complex vector field and polarization field can condense in low tempera- 
tures in an AdS RN black hole background respectively, this model gives a wide possibility 
to investigate the influence between p-wave superconductivity and spontaneous magne- 
tization. We take 7’.9 and Teo as the critical temperatures of py and h, when a = 0. 
Depending on the values of them, the p-wave superconudcting order or ferromagnetism 
will appear first. The interesting question is whether the other phase transition can still 
happen. 


7.3 Coexistence of superconductivity and ferromagnetism 


The first case we will consider is Teo > Teco, i.e., the ferromagnetic phase appears first. 
The equation for c in equations (152) shows that c 4 0 if h 4 0. So there isn’t a phase 
such that {h < 0, ps # 0, py = 0}. When temperature is decreased to lower than 
Too, five kinds of phases may appear. They are phase A {h = ps = py = 0}, phase B 
{Pz = py = 0,h < 0}, phase C {pz # 0,h = py = 0}, phase Dı {ps = py # 0,h < 0} 
and phase D2 {ps = —py # 0,h < 0}, corresponding to normal phase, pure ferromagnetic 
phase, pure p-wave superconducting phase and two kinds of superconducting ferromagnetic 
phases, respectively. 

In this case, whether the p-wave superconductivity can appear depends on the sign of 
interaction of magnetic moment of the complex vector field, i.e., the sign of O. The possible 
phases and the physical favored phase in different temperature regions are summarized in 
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table 1. In the case of © > 0, there is a critical temperature Tse between Too and T 3:0, 
lower than which, the p-wave superconductivity can appear from ferromagnetic phase and 
the system will show ferromagnetism and superconductivity both. In addition, the criti- 
cal temperature for superconductivity is increased rather than decreased by spontaneous 
magnetization. This promotion is enhanced by increasing of interaction strength between 
complex vector field and antisymmetric tensor field. Numerical results imply that the mag- 
netism and superconductivity can coexist even in the zero temperature limit. However if 
© < 0, the p-wave superconducting state can not appear and the system will only be in a 
pure ferromagnetic state. 


Phases in the case of Teo > Tso 


temperature T > Teo | Tse < T < Too T < Toe 
Possible A A,B A, B, Di, Do, C(if T < Tso) 
Physical(© > 0) A B Dı 
Physical(© < 0) A B 


Table 1: The possible and physical phases in the case of Too > Teco. Phase A is {h = pz 
Py = 0}. Phase B is {h < 0, ps = py = 0}. Phase C is {h = py = 0, ps #0}. Phase D; is 
{h <0, px = py #0}. Phase Də is {h < 0, ps = —py # 0}. Table taken from ref. [154]. 


The other case is Tsco > Too, i.e., the case where the p-wave superconducting phase 
appears first. When Too < T < Tæ, the equations (152) show there may exist three 
kinds of p-wave superconducting phases. One is the usual p-wave superconducting phase 
C({h = py = 0, ps # O}), the other two are new superconducting phases denoted as Fy 
with {h < 0, pr = py # 0} and Fa with {h < 0, pr = —p, # 0}. The magnetization in two 
phases Æ; and E» is induced by the p-wave pair rather than been produced spontaneously, 
which is different from case in phases Dı and Dg. 


Phases in the case of Tsc > Teo and O > 0 


Temperature | T > Tso | Too < T < Tico T < Too 
Possible A A, E, C A, E, C, B 
Physical A Ey 


Phases in the case of Teo > Tco and O < 0 


Temperature | T > Tsc | Too < T < Tico | T < Teo 
Possible A A, C A, C, B 
Physical A C B 


Table 2: The possible and physical phases in the case of Tsco > Too. Phase A is {h = pz 


py = 0}. Phase B is {h < 0, px = py = 0}. Phase C is {h = py = 0, pz #0}. Phase Æ is 
{h <0, pr = py #0}. Phases Ey is {h < 0, pz = —py # 0}. The tables were taken from 
ref. [154]. 


Numerical results show that situations also depend on the sign of ©. All the results 
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are summarized in table 2. If © > 0, with decreasing the temperature, the system will 
transit into phase E,, where p-wave superconductivity with a kind of induced magnetism 
appears. The superconductivity and magnetism appear both, however, it should better 
be called a magnetic superconducting phase rather than a ferromagnetic superconducting 
phase, because the magnetic moment is not spontaneously produced and proportional to 
Tso — T rather than /T-9 — T near the critical temperature (see figure 38). If © < 
0, the system will be in the pure p-wave superconducting phase without magnetism if 
temperature is less than Tsco. When temperature is lower than Teco, the system will transit 
into the pure ferromagnetic phase from the p-wave superconducting phase. Therefore the 
ferromagnetism and superconductivity can not coexist if O < 0. 
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Figure 38: The behaviours of N near the critical temperature in the phases Dj; (left) and 
FE, (right). Here m? = —3/16,m3 = —3,J = —1,0 = 1 and a = 0.1. In the left one, 
q = 1.4. In the right one, q = 1.4. The plots were taken from ref. [154]. 


Tables 1 and 2 show that the ground state near zero temperature limit only depends 
on the sign of ©. These can be understood in a physical manner if we pay more atten- 
tion to this phenomenological parameter in (150), where it was introduced to describe 
the self-interaction between the magnetic moments of complex vector field. The case of 
© > 0 means that the p-wave pair will attract the one with the same magnetic moment 
direction and repulse the one with the opposite magnetic moment. Under the influence of 
spontaneous magnetization, the magnetic moment of p-wave pair will tend to align along 
the direction of spontaneous magnetization. As a result, p-wave pair and spontaneous 
magnetization would be enhanced by each other and therefore survive. However, if © < 0, 
the p-wave pair will repulse the one with the same magnetic moment direction. So the 
p-wave pair will align without net magnetism and the system is in a pure p-wave super- 
conducting phase in the region where superconductivity dominates. When T < Too, the 
ferromagnetism can appear, which tends to make p-wave pairs have same direction. But 
the p-wave pairs with same magnetic moment direction will repulse each other, which leads 
that the p-wave pair is not stable and will be de-paired. So the system can only be in the 
ferromagnetic phase. 
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8 Conclusion and Discussion 


Due to the strong/weak duality characteristic of the holographic correspondence, it pro- 
vides us with a powerful tool to study the properties of strongly interacting systems by a 
weakly coupled gravity theory with one extra spatial dimension. Although the underlying 
dynamics which govern the dual field theory and the gravity are apparently different, as we 
have shown, in the framework of holography quantum computations in the dual (strongly 
coupled) field theory can be translated into classical calculations in the bulk, where one 
can just solve differential equations with suitable boundary conditions. Within this frame- 
work, holographic correspondence is considered as a hopeful approach to understand the 
properties of strongly correlated electron systems. 

The bulk gravitational models that we have reviewed are some phenomenological mod- 
els. In such bottom-up approach, the gravity duals were constructed using the minimal 
set of fields that captured the essential dynamics. They just involve gravity interacting 
with an effective U(1) gauge field and a charged field serving as the order parameter. We 
have a lot of the degrees of freedom to choose the form of interactions as well as the value 
of couplings.’ Nevertheless, one has seen that those simple models would describe dual 
superconductors rather well. Some interesting features have been uncovered. Let us take 
holographic p-wave models as an example. For the SU(2) Yang-Mills model (42), the con- 
ductivities are strongly anisotropic in a manner which is suggestive of a gap with nodes. 
The low-lying excitations of the normal state have a relaxation time growing rapidly as the 
temperature is lowered, which agrees with the absence of impurity scattering. For the sec- 
ond model (57), it has been found that the vector condensate can be induced by an applied 
magnetic field, and the condensation of the charged vector operator forms a vortex lattice 
structure in the spatial directions perpendicular to the magnetic field. Going beyond the 
probe approximation, the model displays a rich phase structure. In terms of temperature 
and chemical potential, the complete phase diagrams have been constructed for the con- 
ducting phase, insulating phase and their corresponding superconducting phases and some 
new phase boundaries are revealed. The Maxwell-vector model is a generalization of the 
SU(2) model in the sense that the vector field has a general mass and gyromagnetic ratio. 
The third model (76) realizes a p-wave superconducting phase by involving a charged two- 
form in the bulk. The p-wave states exhibit a helical structure and some of them display 
the phenomenon of pitch inversion as the temperature is decreased. The ground state of 
the condensed phase has zero entropy density and exhibits an emergent scaling symmetry 
in the IR. 

It is clear that the key ingredient in constructing a gravitational dual of a superconduc- 
tor is to find an instability which breaks a U(1) symmetry, e.g., at low temperatures and 
causes a condensate to form spontaneously. One may ask whether those phenomenological 
bulk duals of superconductors are just a Ginzburg-Landau description. The answer is ex- 
actly no. Let us stress two key differences. First, the instability in the Ginzburg-Landau 


28In principle, the arbitrary can be fixed by embedding the bulk model into some low energy effective 
theory of string/M theory. 
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model must be put in by hand, while it arises naturally in holographic setup. Second, the 
Ginzburg-Landau model is only valid near the transition point, whereas the gravitational 
description can characterize the whole dynamics. For a given bulk action, scanning through 
values of model parameters corresponds to scanning through many different dual field the- 
ories. In that sense, a simple holographic model has a kind of universality, i.e., the results 
may be true for a large class of dual field theories, quite insensitive to the details of their 
dynamics. Another confusion is that we realized the spontaneous breaking of a continuous 
U(1) symmetry in (2+ 1) dimensions at finite temperature, in apparent contradiction to 
the Coleman-Mermin-Wagner theorem. The cure is that the large N limit evades the the- 
orem as fluctuations are suppressed. It would be interesting to discuss the effect of bulk 
quantum corrections which correspond to 1/N corrections in the dual field theory [157]. 
Finally, although the hair breaks a local U(1) symmetry in the bulk, according to the 
dictionary, the dual system consists of a condensate breaking a global U(1) symmetry. 
On the other hand, the onset of superconductivity is characterized by the condensation of 
a composite charged operator spontaneously breaking U(1) gauge symmetry. So strictly 
speaking, what one has realized is a dual theory of superfluid [158, 159, 160, 161] rather 
than superconductor. However, in the limit that the U(1) symmetry is “weakly gauged” 
one can still view the dual theory describing a superconductor. 7° 

Throughout this brief summary we have been mainly concerned with static and ho- 
mogeneous case and focused on some basic aspects. This is a rapidly devolving field, 
due to the limitation of length, we are not able to give more details for many interest- 
ing developments, such as introduction of momentum dissipation (to break translational 
symmetry) [164, 165, 166, 167, 168, 169, 170, 171, 172], construction of holographic Joseph- 
son Junction [173, 174, 175, 176, 177, 178, 179, 180, 181, 182], and investigation on dy- 
namics for far-from equilibrium state [183, 184, 185, 186, 187, 188, 189, 190, 191, 192]. 
The analysis for the most part has been done numerically. To explore the properties of 
holographic superconductors using analytical techniques can be found, for example, in 
refs. [193, 194, 195, 196, 197, 198, 199, 200, 201, 202, 203, 204]. Optimistically the growing 
literature based on holographic duality might shed some light on the understanding of 
mysterious phenomena and eventually microscopic origins of strongly correlated supercon- 
ductivity. 

The applications of the holographic correspondence are still going on. It was writ- 
ten by G. T. Horowitz and J. Polchinski [205] that we find it difficult to believe that 
nature does not make use of it, but the precise way in which it does so remains to be 
discovered. In addition to holographic superconductors, the holographic approach has 
been used to understand some other aspects of condensed matter physics, including (non- 
)Fermi liquids [206, 207, 208, 209, 210], quantum Hall effect [211, 212, 213], strange 
metals [214, 215, 216, 217], topological insulators [218, 219, 220], Hubbard model [221] 
and so on. A major application using holographic duality is to describe quantum chro- 


2°In fact, most of the condensed matter theories do not include dynamical photons, as their effects are 
usually small. For example, in the BCS theory electromagnetic field is often introduced as an external 
field. The possibility of introducing dynamical gauge fields in holographic superconductors was discussed 
in refs. [162, 163]. 
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modynamics (QCD), especially for the quark gluon plasma produced in particle accel- 
erators. It is referred as AdS/QCD or holographic QCD, which has been widely stud- 
ied [222, 223, 224, 225, 226, 227, 228, 229, 230, 231, 232, 233, 234, 235]. Another emerging 
subject is the fluid/gravity correspondence, which translates problems in fluid dynamics 
into problems in general relativity [236, 237, 238]. Readers who are interested in those 
exciting achievements are encouraged to consult those relevant references. 
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